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MASSACHUSETTS 


ABSTRACT 


Two  different  procedures  for  analytically  modeling  the  effects  of  the  Moon’s 
direct  gravitational  force  on  artificial  Earth  satellites  are  discussed  from  theoretical 
and  numerical  viewpoints.  One  is  developed  using  classical  series  expansions  of 
inclination  and  eccentricity  for  both  the  satellite  and  the  Moon,  and  the  other 
employs  the  method  of  averaging. 

Both  solutions  are  seen  to  have  advantages,  but  it  is  shown  that  while  the  former 
is  more  accurate  in  special  situations,  the  latter  is  quicker  and  more  practical  for 
the  general  orbit  determination  problem  where  observed  data  are  used  to  correct 
the  orbit  in  near  real  time. 
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1.  INTRODUCTION 


Analytical  theories  describing  the  perturbations  of  artificial  satellites  influenced  by  the  direct 
gravitational  force  of  the  Moon  typically  involve  many  simplifying  assumptions,  for  which  there 
is  little  comment  in  the  literature  as  to  their  effect  upon  the  overall  accuracy  of  the  resulting 
theory.  For  instance,  the  potential  which  gives  rise  to  the  disturbing  function  is  fundamentally 
of  two  degrees  of  freedom  involving  a  complicated  coupling  of  the  motion  of  the  Moon  with  that 
of  the  satellite.  Since  the  disturbing  function  contains  an  exact  generating  function  for  Legendre 
polynomials,  the  problem  is  simplified  somewhat  by  expanding  in  terms  of  these  polynomials  and 
truncating  to  include  just  the  second,  or  sometimes  also  the  third,  harmonic  (where  the  first  is 
absent).  In  order  to  simplify  the  problem  further,  the  motion  of  the  Moon  is  assumed  to  be  a  pure 
Kepler  mean  elliptical  orbit  with  constant  linear  rates  for  the  angular  quantities.  Of  course,  this 
differs  substantial!}  from  the  true  motion  for  the  Moon.  Finally,  in  order  to  perform  a  completeh 
analytical  integration  of  the  equations  of  motion  for  the  satellite,  further  simplifying  assumptions  are 
imposed,  such  as  ignoring  some  of  the  less  dominant  coordinates  over  the  integration  or  truncating 
the  disturbing  function  to  simplify  the  solution. 

In  this  work,  two  philosophically  different  solutions  will  be  presented  and  compared  in  both 
a  theoretical  and  numerical  framework.  It  is  the  intent  for  the  numerical  examples  to  highlight 
the  simplifying  issues  mentioned  above.  The  first  method  to  be  presented  will  be  called  the  "series 
expansion  method"  and  follows  the  development  of  Kaula  1  and  Giacaglia  [2j.  We  will  refer  to 
this  solution  as  the  (SEM)  solution.  The  second  solution  is  developed  using  a  method  of  averaging 
and  follows  the  work  of  Kozai  3  .  Sridharan  and  Seniw  [4],  and  Hujsak  [5j.  We  will  refer  to  this 
solution  as  the  (MAY)  solution.  The  third  harmonic  is  included  in  the  (MAV)  solution  discussed 
in  this  report  (which  is  given  in  its  entirety  in  the  appendix)  since  it  tends  to  be  significant  in 
accurately  modeling  lunar  perturbations  of  high-altitude  satellite  orbits. 

The  gravitational  force  from  a  third  body  acting  on  a  point  source  can  be  described  (assuming 
a  reference  system  defined  by  the  center  of  the  Earth)  using  Newton’s  law  of  gravitation  by  the 
potential 

_  Gmi  (7/70. (ri.  •  r) 

Rk  = - - - -f — l.i 

pk  n 

where  pk  —  r  —  r*  is  the  magnitude  of  the  vector  difference  of  the  vector  r  pointing  to  the  object 
and  the  vector  r^.  pointing  to  the  third  body,  Girtk  is  Newton’s  gravitational  constant  times  the 
mass  of  the  third  bodv.  and  r*  denotes  the  magnitude  of  r*..  For  the  Moon,  the  size  of  Giiik  can 
be  approximated  by 

Cm,  =  N2kal 

where  ^  1.59  x  10“  Drev2/day"  and  Qk  is  the  semimajor  axis  of  the  lunar  orbit  about  the  Earth. 
Using  the  generating  function  for  Legendre  polynomials  in  1  / pk,  we  arrive  at  the  disturbing  function 
for  the  Moon,  which  can  be  written  as 
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(1.2) 


Rk  = 


a- 


E(-)'^(coS  Vk), 

*  —  *  r  i 


where  r  is  the  range  to  the  object  and  v1*  is  the  geocentric  elongation  of  the  satellite  from  the 
Moon  (i.e..  rr^cosvk  =  r  r^.  ).  It  is  worth  noting  that  the  effects  of  the  Moon  are  as  large  or 
larger  than  the  flatness  of  the  Earth  perturbation  for  high-altitude  orbits,  but  are  much  smaller  for 
low-altitude  satellites. 

Kozai  [3],  Sridharan  and  Seniw  [4],  and  Hujsak  [5j  all  truncated  the  disturbing  force  to  /  =  2. 
but  the  development  of  Kaula  1  and  Giacaglia  [2^  is  shown  for  all  /  >  2.  Also  Musen,  Bailie, 
and  Upton  [6]  included  the  third  harmonic  in  their  analysis;  however,  the  ephemeris  of  the  Moon 
in  their  study  relies  on  a  Kepler  element  set  referred  to  the  equator  of  the  Earth.  This  was  also 
done  in  Kaula's  treatment  1  .  It  was  noted  by  Kozai  ^3,7  that  it  is  more  desirable  to  develop  a 
theory  which  relies  on  a  Kepler  ephemeris  for  the  Moon  referred  to  the  ecliptic  plane  since  then  the 
inclination  is  roughlv  constant  and  the  longitude  of  the  right  ascending  node  can  be  approximated 
by  a  simple  linear  function  of  time.  The  work  presented  in  this  report  assumes  this  reference  system 
for  the  lunar  ephemeris  while  the  satellite  elements  are  referred  to  the  equator.  This  will  involve 
some  kind  of  geometrical  rotation  in  both  theories. 
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2.  THE  SERIES  EXPANSION  SOLUTION 


The  (SEM)  solution  has  the  primary  advantage  of  allowing  a  simultaneous  integration  of  the 
two  fundamental  frequencies  involved  in  the  problem,  the  mean  motion  of  the  satellite  and  the 
mean  motion  of  the  Moon.  As  an  added  benefit  of  the  following  analysis,  the  secular  rates  of  the 
argument  of  perigee  and  right  ascending  node  of  the  satellite  and  the  Moon  can  also  be  included  in 
the  integration  at  the  same  time.  The  idea  is  classical  in  principle  and  transforms  the  disturbing 
function  so  that  the  faster  angular  variables  are  converted  to  explicit  functions  of  time  on  the 
right-hand  side  of  the  system  of  differential  equations  so  that  a  direct  solution  is  possible.  The 
development  of  the  solution  is  complicated  and,  for  the  most  part,  quite  handily  displayed  in  2  . 
However,  there  are  some  algebraic  errors  at  crucial  junctures  there,  and  so  a  complete  development 
will  be  presented  here  with  some  of  the  more  evident  mathematical  details  omitted. 

We  can  express  (1.2)  using  Legendre's  addition  theorem 

A-  =  II 

77) )!  «k  a  rk 

Pfn(smtf)Pin(s\n6)cosm(a  -  o').  (2.1) 

where  a.  <L  o'.  6'  are  the  right  ascension  and  declination,  referred  to  the  equator,  of  the  satellite  and 
the  Moon,  respectively: 

1  if  m  =  0 
”  <  2  if  m  =*  0 


a,  ak  denote  the  semimajor  axes  of  the  satellite  and  the  Moon,  respectively;  and  P™ (.r)  denotes  the 
associated  Legendre  function  of  degree  /  and  order  m.  If  we  temporarily  define 

=  Gm (/  ~  ™ ) ^ 1  ( — )'(-)'( ^ )/+ 1  P;m(sin  S') 
ak(l  4-  ?n)!  ak  a  rk 

and  let  C™  =  A™  cos  77m'  and  S™  =  A™  sin  77m',  then  we  can  write 


Rk  —  y  y  P/m(sin^)[C/”  cos  777 o  -f  S™  sin  ma]. 

/> 2  771  =0 

One  can  then  apply  Kaula’s  inclination  functions  1  and  express 

Pa-  =  X!  Flmp(I) 

l> 2  m=0  p—0 
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_ 
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_ 
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\ 

U7 
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7  —  77?  odd 

L  l  J 
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where  0lmp  =  (l  -  2 p)(u:  +  /)  +  mQ  and  u;,  /.  and  f2  denote  the  argument  of  perigee,  true  anomaly, 
and  right  ascending  node  of  the  satellite,  respectively,  and  Fimp(I)  denotes  a  function  of  the  satellite 
inclination  I . 

In  case  that  Kepler  elements  are  chosen  to  represent  the  Moon’s  ephemeris  and  that  they  are 
referred  to  the  ecliptic  plane  (rather  than  the  equator  as  in  the  satellite  case),  a  rotation  of  the 
spherical  harmonics  for  the  Moon  is  needed  [2,  Section  3  .  We  can  write 

i 

Pj71  (sin  fif)  exp  mm']  —  A™'s Pf  (sin  6k)  exp^sa*], 

s=—l 


where  A™’5  is  a  function  of  the  ecliptic  inclination  t,  and  <5*.,  a*.  denote  the  declination  and  right 
ascension  of  the  Moon  referred  to  the  ecliptic  plane,  respectively.  The  expression  for  A™'s  below  is 
taken  directly  from  [8], 


\  U  ^)*  r-/  \  ^  rrrn.St  \ 

A1  (/)  =  exp^m  ■  s)  2 


where 


min  l—s.l—m 

ur(>)  =  (-i  r-  e 

r-maxjO.- (m-rs) 

J  cosm+s+2r(^)sin-rn-^2(,-r)(^). 

Noting  that 


/  —  m 
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p,-‘  (*)-(■ 
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we  can  define 
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and 
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so  that 


cr+isr  =  (i)mY.A7'spi(^h) 


s= 0 

{C,mscos[s(a*  -  |)]  +  iSjn,s  sin[s(aA.  -  ^)]}. 

m  _  _  .... „  { m  —  1 ) 

Since  if  m  is  even  (i)m  =  (  —  1)“  and  if  m  is  odd  ( i)m  =  i(  —  1)  2  ,  it  follows  that 


(  i\k\ST  \-m,sDsi  ■  c  ^m,.s  r  ,  *  ,y  f  CJ”  if  77?  is  ev 

(-1  )‘Z-A  Pf(sm6k)C,  cos[5(qi.  - -)1  =  \  1  . 

JTJ  2  y  5™  if  77 ?  is  od 


en 

odd  : 


and 


(-I)*1  H  A?  ' Pf  (sin  tk)SP's  sur s(ak.  -  ?-)} 

s  —  0  1 


Sfl  if  m  is  even 
- C jm  if  m  is  odd 


where  k\  =  y  .  the  greatest  integer  part  of  y.  Expanding  further  in  terms  of  Kaula’s  inclination 
functions  once  again  (this  time  in  terms  of  the  inclination  Ik  of  the  Moon  referred  to  the  ecliptic 
plane),  we  can  write 


and 


/  / 


(-i)*’ZE-4  ?'sr,sq(h) 


s~Q q=0 


c; 


0 

-yW,S 


"i  l—s  even 


sin  9\ 


Isq 


J  l—s  odd 


s^ttx  ,s 
L  / 
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COS  6\ 


Isq 


l  —  s  odd 


■i 


C;m  if  m  is  even 
S™  if  m  is  odd 


/  / 


(-i 


5=0  <7=0 
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r»m,5 

* l 


-il-s  even 


-sp"' 


l—s  even 


COS  0j 


Isq 


l—s  odd 


sin  #is, ]  >  = 


J  l  —  s  odd 


5;m  if  m  is  even 

—  C™  if  77?  is  odd 


(2.2) 


(2.3) 


where  0J  =  U  ~  2 q)(^'k  +  fk)  +  s(^k  ~  f )  and  a;*.  /*,  and  fl*  denote  the  argument  of  perigee, 
the  true  anomaly,  and  the  right  ascending  node  of  the  lunar  orbit  referred  to  the  ecliptic  plane, 
respectively. 

Since  algebraic  errors  appear  in  the  mathematics  leading  to  equation  (15)  in  [2],  equations  (15 
to  19)  there  are  false;  but  the  form  of  (19)  is  desirable,  and  so  the  mathematics  which  follow  will 
be  aimed  toward  the  objective  of  producing  a  similar  form.  Suppose  m  is  even  and  /  —  rr?  is  even, 
then 
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I  l 

^  =  E  E  E  Flmp(I)  (Qm  CO s6lmp  +  S™  sin 0lmp)  , 

l> 2  m=0p=0 


and  with  some  tedious  algebra,  we  arrive  at 

l  i  l  l 

fi*  =  EEEE  E(-1)^r’'^p(/)f’;„(4)etapJ„  (2.4) 

l>2  m=0  p— 0  5=0  q=0 


Qlmpsq  —  2 


i-iyu?'-* cos (6lmp  +  e\sq)  +  UP'3 cos{dlmp  -  0'hq] 

(-1  r0)m’-s  sin(0lmp  +  6'Ug)  -  U™'s  s\n(6imp  -  6'lsq) 

If  m  is  even  and  /  -  m  is  odd.  then  (2.4)  is  true  with  @impSq  expressed  as 

1  (-1  yU?-'*  sin (0lmp  +  O'uq)  +  Ur  sin (6lmp  -  0'hq)  1 eve" 

-(-1  rur~s  cos(elmp  +  e'lsq)  +  urs  cos(elmp  -  e'lsq)  j;_, odd 

If  rn  is  odd  and  /  —  m  is  even,  then  (2.4)  is  true  with  &impsq  expressed  as 

- 1  )*</”’- sin(»,„p  +  «;„)  +  (7,"'*  sin(«tap  -  »;„)  1 

-iyif- cob(0,mi,  +  «;„)  +  vr'c 05(9,™, -  d;„)  J,., m 


Qlmpsq  ~  2 


Qlmpsq  —  2 
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If  m  is  odd  and  l  —  m  is  odd.  then 


^hmpsq  —  2 


(2.4)  is  true  with  @impsq 

(-l)‘t/,"'-cos(#,mp  +  9;,J)-£/™’'cos(e, l'""v" 


Employing  the  identity  cos(x  —  f )  =  sinx.  we  can  eliminate  the  need  for  the  case  decisions 
on  the  indexing  parameters  by  using  the  following  scheme.  Let  ys  =  §  —  |||||  so  that  ys  —  0  if  s  is 
even  and  ys  —  ^  if  s  is  odd.  By  subtracting  ys7T  from  all  of  the  angular  arguments  of  the  sines  and 
cosines  in  the  above  expressions,  we  can  write 

e,mpsq  =  \{(-l)k2U™'-scos(Olmp  +  0'lsg-ysTr) 

+{~l)k3U™'*  cos(6lmp  -  e'lsq  -  ys:r)}  .  (2.5) 

where  one  solution  for  A'2  and  can  be  expressed  as 


A’2  =  t(m  +  s  —  1)  +  1 
A’3  =  t(m  +  s). 
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with  t  =  (/  —  l)(mod  2)  (that  is,  t  =  0  if  Z  —  1  is  even  and  t  =  1  if  Z  —  1  is  odd).  Substituting  the 
proper  definition  for  A™'s  and  the  right  side  of  (2.5)  for  0/mp5(r  we  can  rewrite  (2.4)  as 


Rk 


iiii 


EErZD-i)* 


a  ,  r  ,  at  ,+ j 
2a*.(/  +  m)!  a*  a  r*  ' 


F,mP(I)F,sq(Ik)  {(-l)kiUjn'~s cos(elmp  +  e\sq  -  ysir) 

+(-l)kiU™'s cos(6hnp  -  8{sq  -  ysir ) }  . 


(2.6) 


Now  Giacaglia  proceeded  to  average  (2.G)  over  the  mean  anomaly  of  the  satellite  to  produce  the 
secular  and  long  periodic  parts  of  the  disturbing  function  together  with  the  short  period  corrections, 
and  a  solution  is  certainly  possible  by  this  route,  but  since  the  real  advantage  of  (2.6)  lies  in  the 
fact  that  it  paves  the  way  for  simultaneous  integration  of  the  angular  rates  for  the  satellite  and 
the  Moon,  we  will  defer  averaging  techniques  to  the  next  section  where  all  averaging  is  performed 
on  the  outset  directly  from  (1.2).  Therefore,  we  will  skip  to  (117)  in  2],  where  the  true  anomalies 
of  the  satellite  and  the  Moon  are  converted  to  their  respective  mean  anomalies  via  the  Hansen 
coefficients,  which  satisfv  the  equation 

(-)n  exp(imf)  =  jr  Hn™r(e)exp[i{m  +  r)M}. 

a  — 


The  final  form  of  the  series  expansion  of  R *  then  becomes 

l  l  l  l  oc  oc 


Rk 


■  -  ZZZ'L'E  Z 


l>  2  771=0  5=0  p=0  <7  =  0  j- -oc  r=—oc 


F,mp(J)Flsq(h)H-_ 


rU-2P 


{(-1  )k2U?''~  ‘cos (8,mpj  +  e'lsqr  -  ysn ) 

+(-l  )k3U™-s  cos (d,mpj  -  e'lsqr  -  j/stt)}  .  (2.7) 

where  6,mpj  =  (/  -  2p)x  -r  (/  -  2p  +  j)M  +  mfi  and  8',sqr  =  (/  -  2q)xk  +  (/  -  2q  +  r)Mk  +  s(Qk  -  § ). 

The  Lagrange  equations  of  motion  of  the  six  Kepler  elements  are  listed  in  [9,  p.  29  and  are 
reproduced  below  for  the  convenience  of  the  reader  : 


2  dR 
nadM 

(1  -  c2)  OR  _  \_J^£dR 
na2e  0M  na^e  dx 


—  cos  I 
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dR  _  \  T^dR 
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na2\  1  -  e2  sin  I 


COS 


dR' 

dft 


_ _ _1 _ dR 

na2\  1  —  e2  sin  I  91 


(1  -e2)9R  _  2_9R 
na2e  de  na  da 


If  the  mean  motion  n  is  approximated  as  n  =  r?o(l  —  where  6a  is  the  perturbation  of  the 

semimajor  axis  due  to  a  above,  and  R k  is  substituted  for  R  in  the  Lagrange  equations,  then  a 
complete  integration  of  the  system  is  possible  provided  a,  a*,  e,  e /,  /*,  t,  .A/,  a?*,  Q,  fl*  are  all 

held  fixed  and  constant  throughout  the  integration.  Numerical  tests  show  that  this  approximation 
is  improved  over  the  long  term  if  the  mean  values  for  these  constants  are  used  instead  of  the  initial 
value  conditions.  Since  this  type  of  solution  is  explained  in  detail  in  [2],  the  mathematical  details 
of  this  straightforward  integration  will  not  be  written  out  here.  All  of  the  work  in  this  section 
carries  over  routinely  for  describing  the  motion  disturbed  by  the  Sun's  gravitational  force,  but  the 
expansion  is  simpler  to  express  since  the  rotation  of  the  spherical  harmonics  is  not  needed  (because 
the  Kepler  elements  for  the  solar  ephemeris  can  remain  in  reference  to  the  equatorial  frame). 


One  final  note  should  be  made  in  this  section  regarding  the  (SEM)  solution.  The  expansion 
involving  the  Hansen  coefficients  is  valid  for  all  eccentricities  less  than  one,  although  the  convergence 
of  these  infinite  series  is  slow  for  high  eccentricity.  The  interested  reader  is  encouraged  to  consult 
[10~  for  a  detailed  study  of  this  point.  Moreover,  it  is  well  known  that  the  sum  over  j  cannot  be 
truncated  to  just  a  few  values  if  the  eccentricity  is  high,  although  the  series  over  r  can  be  truncated 
to  —3  <  r  <  3  involving  the  eccentricity  of  the  Moon  (e*  ~  .054).  These  issues  make  the  (SEM) 
solution  extremely  impractical  for  use  with  all  satellites.  Even  if  a  satellite  has  small  eccentricity, 
a  truncation  of  Rk  in  (2.7)  to  /  =  2.3  with  —  3  <  j  <  3  can  involve  literally  thousands  of  terms  in 
any  one  solution.  This  has  the  sobering  effect  of  causing  the  analytical  solution  outlined  here  to  be 
actually  slower  than  a  full-blown  numerical  integration  over  (not  too)  small  time  scales,  say,  less 
than  five  days.  The  argument  here  is  a  strong  motivation  for  developing  a  solution  by  a  method 
which  does  not  resort  to  any  expansions  save  the  first  one  involving  Legendre  polynomials,  and 
which  requires  the  same  computation  power  no  matter  the  satellite  parameters.  Nevertheless,  the 
(SEM)  solution  is  mathematically  pleasing  since  it  displays  the  entire  solution  for  all  harmonics 
in  a  single  compact  form,  and  the  dissection  of  the  perturbations  into  secular,  long  periodic,  and 
short  periodic  parts  can  be  readih  made  and  studied. 
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3.  THE  METHOD  OF  AVERAGING  SOLUTION 


The  (MAY)  solution  to  be  described  in  this  report  is  derived  by  an  entirely  different  approach 
than  the  (SEM)  solution.  For  starters,  averaging  techniques  look  for  closed  form  expressions  of  a 
solution  without  resorting  to  infinite  series  expansions  of  the  type  outlined  in  the  previous  section. 
In  order  to  meet  this  objective,  the  differential  equations  are  transformed  so  that  the  faster  periodic 
variables  are  eliminated  and  the  resulting  equations  become  simpler  to  integrate.  Averaging  meth¬ 
ods  are  particularly  useful  in  situations  where  more  than  one  degree  of  freedom  is  present,  and  in 
the  case  at  hand,  it  is  assumed  that  the  mean  anomaly  of  the  satellite  and  the  mean  anomaly  of  the 
Moon  are  fundamental  periodic  quantities  to  be  averaged.  The  averaging  allows  integration  of  the 
equations  of  motion  by  one  variable  at  a  time,  and  hence,  series  expansions  in  terms  of  eccentricity 
are  avoided.  It  will  be  seen  that  the  (MAY)  solution  produces  useful  results  and  is  appropriate 
for  a  wide  range  of  satellite  applications  because  it  does  not  have  this  inherent  difficulty  with  high 
eccentricity. 

A  few  words  should  be  mentioned  concerning  the  averaging  method  used  here.  In  some  ways, 
the  technique  is  similar  to  averaging  methods  used  in  other  satellite  orbit  problems  (such  as  the 
main  satellite  problem),  but  some  shortcuts  have  been  taken  (when  sufficient  evidence  is  warranted) 
because  of  the  complicated  nature  of  our  problem  at  hand.  In  the  von-Zeipel  method,  for  example, 
the  theory  of  canonical  transformations  is  exploited  and  a  generating  function  is  used  to  transform 
the  old  Hamiltonian  and  old  coordinates  to  a  new  Hamiltonian  with  new  coordinates  such  that 
a  high  frequency  variable  is  eliminated  and  the  differential  equations  become  simpler  to  solve. 
Successive  transformations  lead  to  constant  action  \ariables  and  to  angle  variables  wdth  constant 
linear  rates,  which  become  the  foundation  for  a  complete  solution.  As  long  as  no  resonances 
exist,  this  technique  has  been  used  for  varied  problems  in  mechanics  with  favorable  results.  The 
generating  functions  are  typically  constructed  by  using  Taylor  series  expansions  and  integrating  the 
appropriate  partial  derivative  components,  and  truncation  allows  the  expressions  to  be  simplified  to 
desired  orders.  Generating  functions  define  the  transformation  between  new  and  old  coordinates, 
and  this  process  of  transformation  is  usually  referred  to  by  saying  that  the  old  coordinates  are 
obtained  by  adding  periodic  corrections  to  the  new  variables.  An  excellent  example  of  this  method 
used  in  practice  is  Brouwer’s  solution  of  the  main  satellite  problem  [13],  The  generalized  method  of 
averaging  relies  on  asymptotic  expansions  to  accomplish  the  same  goal,  only  the  canonical  form  is 
not  necessarily  preserved  after  transformation.  It  is  shown  in  [11  how  the  von-Zeipel  method  is  a 
special  case  of  the  generalized  method  of  averaging  (at  least  through  second  order),  and  it  is  argued 
that  many  techniques  which  require  specific  transformations  can  fall  under  the  more  powerful  and 
general  category  of  asymptotic  expansions. 

In  our  use  of  the  method  of  averaging,  a  solution  is  desired  which  is  accurate  to  first  order  only, 
and  so  the  disturbing  function  is  averaged  directlv  in  the  spirit  of  [3].  The  differential  equations 
are  separated  into  their  secular  (those  which  have  the  mean  anomalies  of  the  satellite  and  the 
Moon  absent),  long  periodic  (those  which  have  an  explicit  dependence  on  the  mean  anomaly  of  the 
Moon),  and  short  periodic  (those  wffiich  have  an  explicit  dependence  on  the  mean  anomaly  of  the 
satellite)  parts.  This  is  entirely  consistent  wTith  the  generalized  method  of  averaging  (to  fir^t  order): 
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however,  after  two  transformations  by  averaging  these  anomalies,  it  turns  out  that  the  transformed 
action  variables  are  not  yet  constant.  The  rates  of  the  eccentricity  and  inclination  are  not  zero 
since  the  argument  of  perigee  and  right  ascending  node  of  the  satellite  and  the  Moon  have  not 
been  averaged  out.  Nevertheless,  the  integration  is  still  proceeded  through  in  the  same  spirit  of 
the  generalized  method  of  averaging,  because  the  error  in  so  doing  is  assumed  small.  It  will  be 
shown  numerically  in  the  next  section  that  this  approach  yields  a  tractable  solution  which  is  a  good 
first-order  approximation  for  a  varied  class  of  satellite  orbits. 

The  theoretical  framework  for  the  (MAY)  solution  is  given  in  [4],  and  in  this  section  the 
procedure  will  be  outlined  with  some  technical  details  noted.  The  actual  solution  produced  in  the 
appendix  of  [4]  is  given  for  just  the  I  =  2  harmonic  in  (1.2)  while  the  solution  in  the  appendix 
of  this  report  corrects  some  minor  algebraic  errors  and  also  includes  the  l  =  3  harmonic.  This 
has  proved  necessary  for  high-altitude  satellites  such  as  those  in  synchronous  or  half-synchronous 
orbits.  The  (MAY)  solution  involves  an  enormous  amount  of  algebra,  especially  in  averaging  over 
the  mean  anomaly  of  the  Moon.  This  is  a  mathematical  disadvantage,  since  there  is  no  way  to  rely 
on  the  work  which  produced  the  l  =  2  solution  for  the  /  =  3  solution  as  in  the  (SEM)  case.  Because 
of  the  bulky  amount  of  algebra  involved  in  the  (MAY)  solution,  computer  algebra  is  very  helpful  in 
performing  the  integration  since  it  can  produce  fast,  error-free  computations  with  FORTRAN  code 
also  available.  The  solution  presented  in  the  appendix  was  derived  by  MIT’s  MACSYMA  symbolic 
manipulation  program  and  a  FORTRAN  implementation  on  magnetic  tape  of  either  the  (SEM) 
solution  or  the  (MAY)  solution  is  available  from  the  author  upon  request. 

If  z  represents  a  vector  containing  the  six  Lagrange  differential  equations  describing  the  motion 
of  the  vector  of  satellite  elements  z  under  the  disturbing  force  of  the  Moon  (1.2),  the  (MAY)  solution 
is  started  with  an  average  of  z  over  the  mean  anomaly  of  the  satellite 

1  f2n 

<Z>  =  IT  /  zdM.  (3.1) 

27T  Jo 

In  the  case  at  hand,  it  is  convenient  to  use  the  eccentric  anomaly  E  to  compute  the  indefinite 
integral 


h  =  I  z(£)(l  ~  ecosEjdE 

so  that  the  short  period  corrections  appear  as  a  by-product 
nb  zs  =  /]  —  E{  z). 

where  n  denotes  the  mean  motion  of  the  satellite.  Since  this  is  an  initial  value  problem,  we  can 
remove  the  constant  part  of  this  quantity.  A  second  averaging  is  now  performed  on  (z)  with  respect 
to  the  mean  anomaly  of  the  Moon 


1  f2ir 

«*»-sZ  < 


(z)dMk- 
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and  this  is  calculated  by  converting  (z)  in  terms  of  the  true  anomaly  of  the  Moon  and  computing 
the  indefinite  integral 


h  =  J 


The  long  period  corrections  are  readily  available  as 


nktzi  =  I2  -  A-«z». 


where  n*.  denotes  the  mean  motion  of  the  Moon  Thus  to  update  given  mean  elements  zq  (at  a 
given  epoch  to)  to  a  desired  time  t ,  one  propagates  the  mean  elements 

Z)  =  z0  +  ((£}>(*  -  to)- 

adds  the  long  period  corrections 

Z'2  =  Zi  -r  6z/(Zj  ). 

and  then  adds  the  short  period  corrections 

Z  (/)  =  Z2  +  (3zvS(z2). 

Since  the  equation  for  the  mean  anomaly  rate  M  has  the  mean  motion  n  included  and  (a)  =  0.  it 
is  necessary  in  a  first-order  development  to  include 


3  n6as 
2  ~ 


in  the  equation  of  M  before  any  averaging  of  this  equation  is  performed. 


Some  technical  details  involved  in  the  averaging  mentioned  above  will  now  be  presented.  Let 
(i.j,k)  denote  a  fixed  Earth  centered  inertial  coordinate  system  referred  to  the  equator  (with  i 
lying  in  the  ecliptic  plane  pointing  to  a  vernal  equinox)  and  let  Zq  be  a  vector  of  the  six  initial 
Kepler  elements  referred  to  this  system.  Now.  let  P  denote  the  unit  vector  pointing  to  perigee  of 
the  elliptical  satellite  orbit  and  let  Q  =  denote  the  unit  vector  orthogonal  to  P  in  the  satellite 


plane.  To  complete  a  right-hand  system,  let  u\-  denote  the  unit  normal  to  the  satellite  plane. 
u,v  =  P  x  Q.  It  is  well  known  that 


P  =  ( cos  17 cos uj  —  cos  /sinf2siiiwL',)i  T 

(sin  Q  cos  u;  +  cos  /  cos  Q  sin  u;)j  -f 
(sin  /  sino?)k 


(3.2) 


and 
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u\  =  (sin  ft  sin  7)i  -  (cosftsin/)j  4  (cos7)k< 


(3.3) 


The  unit  vector  lr  pointing  to  the  satellite  at  any  time  can  be  written 
lr  =  cos(/)P  -i-  sin(/)Q. 

and  if  1*.  denotes  a  unit  vector  pointing  to  the  Moon,  for  the  first  averaging  it  suffices  to  use  the 
projections  of  1^.  onto  P.Q.  and  u\  .  These  projections  will  be  denoted  51,  52,  and  53,  respectively, 
so  that  we  can  express 

cos xpk  =  lr  •  1*  =  cos(/)51  4  sin(/)52. 

The  projection  53  =  u.y  1  k  appears  in  the  Lagrange  equations  for  the  rates  of  the  argument  of 
perigee  bb.  right  ascending  node  ft.  and  the  inclination  7  since  =  u^ysinu;  and  cos  7^f  —  = 

-  sin  /|j. 

For  the  second  averaging  over  the  mean  anomaly  of  the  Moon,  some  geometrical  transforma¬ 
tions  are  needed.  From  Kepler  elements  for  the  Moon's  orbit,  we  can  compute  u^.,  the  unit  normal 
of  the  Moon’s  orbital  plane,  and  P* the  vector  pointing  to  the  Moon's  perigee,  as  in  (3.3)  and  (3.2). 
respectively.  If  the  elements  for  the  Moon's  orbit  are  referred  to  the  ecliptic,  then  an  additional 
transformation  must  be  made  to  place  these  vectors  into  the  inertial  reference  frame  (i,j.k).  This 
is  merely  a  plane  rotation  about  the  ecliptic  inclination  t;  for  if  v  =  cii*  4-  c^l  +  03k,  is  any  vector 
referred  to  an  inertial  system  about  the  ecliptic  plane  (with  iL  coinciding  with  i),  then  we  can  write 
v  =  ci i  4-  (C2  cost  -  C3  sin  i )j  4  (co  sin  t  4  C3  cos  t)k  in  the  inertial  reference  frame  about  the  equator. 
As  in  the  case  of  the  satellite  orbit,  if  we  define  Qk  =  u k  x  P^,  then  we  can  express 

1a  =  cos  (A-)  Pa-  +  sin(A)QA-- 

It  is  now  convenient  to  write  the  inclination  of  the  Moon  with  respect  to  the  satellite  plane 


cos  jk  EUy  uk. 


and  compute  the  unit  vector  pointing  in  the  direction  of  the  nodal  line  which  intersects  the  plane 
of  the  satellite  orbit  and  the  plane  of  the  lunar  orbit. 

uk  x  u.v 

u  L  =  — T  • 
sm  jk 

To  complete  a  right-hand  system  for  this  representation  of  the  satellite  plane,  we  define 


u xl  =  u.v  x  uL 


so  that 
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P  =  cos(apk)ui  4-  s\n(apk)usi 


and 


Q  =  —  sm(apk)ui  -I-  cos(apk)uxL> 

where  apk  is  the  argument  of  perigee  of  the  satellite  orbit  with  respect  to  the  Moon's  orbit.  To 
complete  a  right-hand  system  for  this  representation  of  the  lunar  orbit,  we  define 


UAL  =  Ua  x  UL 


so  that 


Pa-  =  -cos(u'k)uL  -  s\n(u'k)ukL 


and 


Qa-  =  sin(  wk)ui  -  cos{u'k)ukL • 

where  wk  is  the  argument  of  perigee  of  the  lunar  orbit  measured  with  respect  to  the  satellite  orbit. 
Note  that  the  nodal  vector  ul  points  to  the  descending  node  of  this  system  for  the  Moon  but  points 
to  the  ascending  node  of  the  system  for  the  satellite. 

It  is  not  difficult  to  use  this  information  to  represent  51,  52,  and  53  in  their  final  form  so  that 
an  averaging  over  the  mean  anomaly  of  the  Moon  can  be  performed: 

51  =  —  cos(apA’)  c.os(u'k  4-  fk)  —  sii \(apk)  cos(jk)  s\n(uxk  4-  fk) 

52  =  sin(apA')  cos(u'k  —  fk)  —  cos(apk)  cos( jk)  s\n(u'k  4-  fk) 

53  =  sin(jA')  sii^u’A'  4-  fk). 

Of  course,  all  of  the  equations  carry  over  verbatim  for  use  with  the  solar  disturbing  function  with 
appropriate  consideration  for  the  values  pertaining  to  the  Sun’s  force  in  place  of  those  for  the 
Moon's  force. 
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4.  NUMERICAL  COMPARISONS 


The  (SEM)  and  (MAY)  solutions  have  been  coded  in  FORTRAN  for  a  Harris  800  series  com¬ 
puter  with  double  precision  (48-bit)  arithmetic.  These  solutions  will  be  compared  numerically  in 
two  settings.  One  is  a  comparison  of  each  analytical  solution  against  a  straight  numerical  integra¬ 
tion  of  the  Lagrange  equations  of  motion,  and  the  second  is  a  test  to  see  how  well  each  solution 
performs  in  an  analytical  orbit  determination  routine  which  fits  (in  the  weighted  least  squares  sense) 
observed  data  to  a  model  of  propagation.  The  observed  data  in  this  case  are  a  set  of  topocentric 
sites  pointing  to  an  object  consisting  of  azimuth,  elevation,  range,  and  range  rate.  The  data  were 
simulated  by  a  special  perturbations  precision  numerical  integration  routine  which  integrates  New¬ 
ton’s  law  of  gravitation  directly  from  (1.1)  (and  therefore  requires  no  truncation  of  a  Legendre 
expansion).  Simulated  data  were  chosen  over  real  data  for  the  following  reasons:  simulated  data 
were  essentially  free  of  noise,  other  forces  can  be  turned  on  or  off  as  desired  so  that  interactions  can 
be  investigated,  and  all  residuals  are  directly  related  to  the  analytical  model  instead  of  possibly 
belonging  to  some  unknown  effect  in  the  data.  An  additional  advantage  to  using  simulated  data 
were  that  various  ephemeris  models  for  the  Moon  can  be  placed  into  the  integration  such  as  the 
simple  mean  Kepler  theory  (which  is  assumed  by  the  analytical  theories),  a  more  accurate  and 
complex  lunar  theory  such  as  the  Hill-Brown  model  12  ,  or  interpolation  from  ephemeris  tapes 
which  have  been  provided  by  the  U.S.  Naval  Observatory.  In  any  case,  all  that  is  needed  in  a 
numerical  integration  of  the  equations  of  motion  is  a  unit  vector  pointing  to  the  Moon  together 
with  its  range  at  any  desired  time. 

All  of  the  tests  below  involve  a  generic  geosynchronous  object  with  parameters:  a  =  42,  1G4.1 
km.  e  —  .001.  I  —  1.5  deg.  ^  =  145  deg.  =  16Gdeg.  and  M  =  25.413 deg  at  the  epoch  day 
280.28046  of  1986.  Considering  the  present-day  satellite  population,  the  effect  of  the  Moon's  grav¬ 
itational  force  is  most  significant  at  synchronous  altitudes,  with  the  obvious  exception  of  super- 
synchronous  objects  with  high  eccentricity  such  as  the  Soviet  Prognoz  series.  Such  objects  reach 
over  half  of  the  distance  to  the  Moon  at  apogee,  and  to  use  this  class  of  satellites  in  numerical 
experiments  with  the  analytical  theories  in  question  would  be  pointless  because  a  truncation  of  the 
force  model  to  just  the  second  and  third  harmonics  is  no  longer  justified.  On  the  other  hand,  a 
geosynchronous  object  is  ideal  for  numerical  comparisons  of  this  nature  because  the  perturbations 
are  significant,  the  eccentricity  is  small  (so  that  a  truncation  of  the  Hansen  coefficients  in  the  (SEM) 
solution  can  be  chosen  at  I  j  <  3).  the  third  harmonic  plays  a  noticeable  role  in  achieving  a  desired 
accuracy,  and  the  class  of  synchronous  objects  is  one  of  the  largest  classes  of  active  satellites  most 
affected  by  the  Moon's  gravitational  force. 

First,  the  (SEM)  and  (MAY)  solutions  were  compared  to  a  fifth-  and  sixth-order  variable  step 
size  Runge-Kutta  numerical  integration.  Recall  that  both  of  the  analytical  solutions  assume  that 
the  ephemeris  for  the  Moon  is  described  by  a  Kepler  ellipse  with  constant  rates  of  change  for  the 
angular  quantities.  In  the  numerical  integration,  it  is  possible  to  insert  any’  ephemeris  model  for  the 
Moon.  \Ye  will  use  two  different  ephemeris  models  for  comparison  in  this  experiment.  The  Kepler 
ephemeris  which  is  used  for  the  analytical  solutions  will  be  referred  to  as  the  “simple  model."  and 
a  more  complicated  and  accurate  model  known  as  the  Hill-Brown  lunar  theory  [12  will  be  referred 
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to  as  the  “Hill-Brown  model.**  The  purpose  in  this  is  that  it  allows  one  to  see  the  effect  in  pure 
propagation  of  simplifying  the  ephemeris  model  so  that  integration  is  possible.  In  each  of  the  plots 
which  follow,  the  effect  upon  the  semimajor  axis  and  mean  longitude  (6 A  =  6u>  +  6M  +  £>Q)  is 
displayed. 

Figure  4-1  displays  the  total  perturbation  in  a  and  A  from  the  Moon's  direct  gravitational 
force.  Both  the  second  and  third  harmonics  were  used  in  the  force  model  and  the  simple  lunar 
ephemeris  model  was  used  in  the  numerical  integration.  Notice  that  the  amplitude  in  the  change  to 
the  semimajor  axis  reaches  2  km  and  that  the  perturbation  of  the  mean  longitude  has  a  strong  sec¬ 
ular  rate.  It  was  observed  (but  not  plotted  here)  that  the  amplitude  of  the  change  in  the  semimajor 
axis  from  the  second  harmonic  in  the  force  model  was  approximately  2  km,  the  amplitude  from  the 
third  harmonic  in  the  force  model  was  approximately  300  m,  and  the  amplitude  from  the  fourth 
harmonic  in  the  force  model  was  approximately  20  m.  The  comparison  of  this  graph  to  the  (SEM) 
solution  is  shown  in  Figure  4-2.  and  the  comparison  to  the  (MAY)  solution  is  shown  in  Figure  4-3. 
It  appears  that  under  these  special  circumstances,  the  (SEM)  solution  exhibits  smaller  error  resid¬ 
uals  by  an  order  of  magnitude.  The  errors  in  the  (SEM)  solution  behave  as  one  would  expect  from 
a  first-order  theory  (about  10“3  to  10“4  of  the  perturbation),  but  the  (MAY)  solution  has  larger 
errors  in  the  short  period  corrections.  These  differences  reflect  upon  the  approximation  to  the  true 
solution  for  each  theon.  however,  they  become  less  pronounced  when  the  Hill-Brown  ephemeris  is 
used  in  the  numerical  integration.  The  error  in  the  (SEM)  solution  for  this  situation  is  displayed  in 
Figure  4-4,  where  one  can  now  see  the  amplitude  of  the  error  in  the  semimajor  axis  reaching  100  m. 


TIME  (In  Days)  FROM  EPOCH 
SEMIMAJOR  AXIS  (km)  ♦  LONGITUDE  {E-1  dag) 

Figure  4-1.  The  lunar  perturbation  of  the  semimajor  axis  and  mean  longitude. 

The  length  of  propagation  is  another  parameter  which  may  tell  more  about  the  effective  ac¬ 
curacy  of  the  two  solutions  In  Figures  4-5  and  4-6,  the  error  in  the  propagation  of  the  semimajor 
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114043  1 


114043-3  114043-2 


Figure  4-2.  The  (SEM)  error  in  the  seiniinajor  axis  and  the  mean  longitude  with  the 
simple  lunar  ephemeris  model  used  in  the  numerical  reference. 


TIME  (In  Days)  FROM  EPOCH 

SEMIMAJOR  AXIS  (m)  ♦  LONQfTUDE  (E-4  dag) 

Figure  4-3.  The  (MAY)  error  in  the  semima jor  axis  and  the  mean  longitude  with  the 
simple  lunar  epheuieris  model  used  in  the  numerical  reference. 
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Figure  4-4  The  (SEM)  error  in  the  semiwajor  axis  and  the  mean  longitude  with  the 
Hill-Brown  lunar  epheineris  model  used  in  the  numerical  reference. 


axis  and  mean  longitude  is  extended  to  100  days  from  the  epoch  for  the  (SEM)  solution  and  the 
(MAY)  solution,  respectively.  This  is  for  comparison  with  the  plots  in  Figures  4-2  and  4-3.  and  so 
the  simple  lunar  ephemeris  model  was  used  in  the  numerical  integration.  Incidentally,  the  integra¬ 
tion  constants,  which  are  different  for  both  theories,  contribute  to  the  error  growth  and  apparent 
biases  displayed  in  these  plots.  This  is  of  no  great  concern  since  when  each  theory  is  fit  to  observed 
data,  these  constants  are  automatically  fit  in  the  process.  These  plots  show  that  even  over  100 
revolutions  of  the  satellite  orbit,  both  theories  are  able  to  predict  the  lunar  perturbations  to  within 
reasonable  tolerances. 

In  many  cases,  an  analytical  theory  is  chosen  so  that  speed  in  computation  can  be  achieved 
for  routine  calculations.  Often,  near  real  time  experiments  are  desired,  and  precision  numerical 
theories  are  not  able  to  meet  this  requirement.  Thus,  it  is  necessary  to  consider  computation  speed 
as  an  important  parameter  in  order  to  determine  the  usefulness  of  an  analytical  theory.  Table 
4-1  summarizes  the  time  (in  milliseconds)  required  for  one  propagation  from  each  solution  on  the 
Harris  H-800  computer,  where  the  force  model  included  just  the  second  harmonic,  or  the  second 
and  the  third  harmonics.  It  is  seen  that  the  (MAY)  solution  is  faster  by  an  order  of  magnitude. 
This  is  not  surprising  since  a  rough  count  of  the  number  of  terms  in  (2.7)  indicates  that  there  are 
3x3x3x3x7x7  =  3,  969  periodic  terms  in  the  expanded  force  model  for  just  the  l  =  2  Legendre 
polynomial,  and  for  the  /  =  2  and  /  =  3  force  model,  there  are  3,969  +  4x4x4x4x7x7  = 
16.513  periodic  terms.  Some  of  these  terms  can  be  neglected  in  a  “smart”  implementation  of  the 
theory  because  they  contribute  little  to  the  overall  expansion,  and  clever  bookkeeping  is  required  to 
optimize  the  running  time.  However,  it  is  clear  that  this  solution  contains  a  lot  of  terms  and  requires 
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TIME  (In  Day*)  FROM  EPOCH 

SEMIMAJOR  AXIS  (m)  +  LONGITUDE  (E-5  deg) 
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Figure  4-5.  The  (SEM)  error  in  the  semimajor  axis  and  the  mean  longitude  with  the 
simple  lunar  ephemeris  model  used  in  the  numerical  reference. 


so  - 


0  20  40  60  80  100 

TIME  (In  Days)  FROM  EPOCH 

SEMIMAJOR  AXIS  (m)  +  LONGITUDE  (E-3  deg) 

Figure  4-6.  The  (MAV)  error  in  the  semimajor  axis  and  the  mean  longitude  with  the 
simple  lunar  ephemeris  model  used  in  the  numerical  reference. 
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the  evaluation  of  hundreds  (and  possibly  thousands)  of  sines  and  cosines.  While  the  (M  AY)  solution 
can  also  have  thousands  of  terms,  there  are  only  a  handful  of  different  trigonometric  functions  to 
evaluate,  and  this  causes  the  overall  running  time  for  this  solution  to  be  much  shorter  than  for  the 
(SEM)  solution  (since  multiplication  and  addition  take  less  time  to  compute  than  the  evaluation 
of  transcendental  functions). 


TABLE  4-1 


Computation  Time  Required  for  One  Propagation 


Force  Model 

(SEM)  solution 

(MAV)  solution 

/  =  2  harmonic  only 

85  ms 

11  ms 

/  =  2  and  l  =  3  harmonics 

256  ms 

27  ms 

The  final  numerical  experiment  involves  fitting  mean  elements  to  numerically  simulated  data 
using  the  (SEM)  solution  or  the  (MAV)  solution  as  a  propagator.  Ten  days  of  observed  data  were 
generated  using  precision  numerical  integration  of  the  special  perturbations  from  the  lunar  force 
model.  The  data  were  propagated  directly  from  (1.1),  and  therefore,  no  expansion  (or  truncation) 
of  Legendre  polynomials  was  employed.  Two  sets  of  data  consisting  of  azimuth,  elevation,  range, 
and  range  rate  from  a  selected  site  location  were  generated,  with  the  first  using  the  simple  ephemeris 
model  for  the  Moon  and  the  second  using  the  more  realistic  Hill-Brown  ephemeris  for  the  Moon. 
The  converged  residuals  are  presented  in  Table  4-2.  It  is  interesting  to  note  that,  in  this  experiment, 
there  is  very  little  difference  in  accuracy  between  the  two  analytical  theories.  It  is  also  clear  that 
a  mean  or  osculating  ephemeris  for  the  Moon  is  not  a  disturbing  factor  in  this  experiment  (as  is 
expected  since  this  is  an  iterative  process).  The  third  row  in  Table  4-2  is  presented  only  to  show 
the  large  residuals  in  range  which  result  when  the  third  harmonic  is  omitted  from  the  force  model 
for  this  particular  object.  The  accuracy  in  range  is  degraded  by  over  300  m. 

In  spite  of  the  results  from  the  direct  propagation  experiment,  the  fitting  experiment  and  the 
time  computation  experiment  demonstrate  that  the  (MAV)  solution  is  more  practical  and  just  as 
accurate  for  most  purposes.  The  assumptions  provide  sufficient  accuracy  to  within  70  to  100  m  in 
range  and  2  mdeg  in  azimuth  and  elevation,  and  it  is  doubtful  that  these  levels  can  be  reduced  by 
any  addition  to  the  theory. 
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TABLE  4-2. 


Converged  rms  Residuals  as  Compared  Against  Ephemeris  Model  in  the  Reference  and 
Force  Model  in  the  Solution 


(MAV)/(SEM) 

Az.  (mdeg) 

El.  (mdeg) 

r  (m) 

dotr  (cm/s) 

Simple  ephemeris  model 
/  =  2,3  harmonics 

1  i 

1 

i  / 1 

70  /  80 

1 

'  1 

Hill-Brown  model 

/  =  2.3  harmonics 

i  / 

1 

i  / 1 

60  /100 

1 

1  1 

—  ■ 

Hill-Brown  model 

/  =  2  harmonic 

1 

2 

2 

i  / 1 

420  /  420 

3 

3 
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5.  SUMMARY 


Two  analytical  solutions  to  the  problem  of  describing  the  direct  lunar  gravitational  effect  upon 
artificial  satellite  orbits  have  been  presented  and  compared  from  both  a  theoretical  and  numerical 
viewpoint.  The  (SEM)  solution  uses  traditional  expansions  about  the  inclination  and  eccentricity 
of  both  the  satellite  and  the  Moon  so  that  a  simultaneous  integration  of  both  motions  is  afforded, 
while  the  (MAY)  solution  employs  the  method  of  averaging  on  the  unexpanded  disturbing  function 
(1.2)  so  that  the  motions  of  the  satellite  and  the  Moon  can  be  integrated  separately.  The  (SEM) 
solution  is  complete  for  all  Legendre  harmonics  in  the  disturbing  function  (1.2),  whereas  the  (MAY) 
solution  is  complicated  to  express  for  more  than  the  second  and  third  harmonics.  The  following 
conclusions  can  be  drawn. 

•  The  (SEM)  solution  is  more  pleasing  mathematically  because  it  provides  a  com¬ 
pact  expression  for  the  complete  problem,  while  the  (MAY)  solution  presents 
a  bulky  package  of  equations  for  one  harmonic  at  a  time. 

•  The  simultaneous  integration  exhibited  by  the  (SEM)  solution  achieves  about 
an  order  of  magnitude  (or  more)  in  accuracy  over  the  (MAY)  solution  w’hen 
comparing  the  perturbation  of  the  Kepler  elements  of  the  satellite  against  a 
numerical  integration  of  the  Lagrange  equations  of  motion  with  an  ephemeris 
for  the  Moon  approximated  by  mean  Kepler  motion.  Under  these  conditions, 
a  ten-day  integration  of  a  geosynchronous  object  exhibits  an  error  of  about 
3  m  and  70  m  for  the  semimajor  axis  from  the  (SEM)  solution  and  the  (MAY) 
solution,  respectively.  The  error  in  mean  longitude  is  approximately  10~5  deg 
for  the  (SEM)  solution  and  10“4  deg  for  the  (MAV)  solution.  The  total  pertur¬ 
bation  in  this  experiment  reached  2  km  for  the  semimajor  axis  and  10~ 1  deg 
for  the  mean  longitude. 

•  The  (SEM)  solution  can  be  more  than  an  order  of  magnitude  slower  in  compu¬ 
tation  speed  for  one  propagation  than  the  (MAY)  solution,  and  if  near  real  time 
execution  is  necessary,  then  the  (SEM)  solution  is  in  serious  jeopardy  of  being 
entirely  impractical  for  this  purpose.  Moreover,  the  execution  time  required 
for  the  (SEM)  solution  increases  significantly  for  objects  wfith  high  eccentricity, 
while  the  (MAY)  solution  shows  no  increase  for  such  objects. 

•  Both  the  (SEM)  and  (MAY)  solutions  perform  the  same  wrhen  asked  to  fit 
simulated  data  for  a  geosynchronous  object  over  a  ten-day  period.  Converged 
residuals  wrere  on  the  order  of  a  millidegree  in  azimuth  and  elevation,  70  m  in 
range,  and  10  min/s  in  range  rate. 

•  The  ephemeris  model  has  little  effect  on  the  converged  residuals  for  both  mod¬ 
els,  and  this  suggests  that,  for  this  application,  the  assumption  of  mean  Kepler 
motion  for  the  Moon  in  the  analytical  theories  is  good  for  a  first-order  theory. 
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•  Both  the  second  and  third  harmonics  in  the  expansion  of  the  disturbing  function 
(1.2)  are  necessary  for  high-altitude  orbits  in  order  to  obtain  an  accuracy  of 
less  than  100  m  in  range. 
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APPENDIX  A 


The  method  of  averaging  solution  for  the  l  =  2  and  /  =  3  harmonics  in  the  disturbing  function 
is  presented  below.  The  solution  was  produced  and  checked  by  MACSYMA,  a  computer  algebra 
manipulation  program  developed  at  MIT.  The  solution  should  be,  for  the  most  part,  free  from 
errors  (although  the  process  of  transferring  the  solution  to  paper  can  involve  human  error).  The 
equations  which  are  reproduced  below  have  been  coded  in  FORTRAN  (also  by  MACSYMA)  and 
numerical  tests  indicate  that  the  the  equations  are  correct.  The  solution  is  written  in  terms  of 
quantities  already  defined  except  for  combinations  of  magnitude  which  will  be  defined  below. 


The  first  averaging  is  over  the  mean  anomaly  of  the  satellite,and  for  the  l  =  2  harmonic  we 
define 


where  is  the  ratio  of  the  mass  of  the  Moon  to  the  mass  of  the  Earth.  For  the  l  =  3  harmonic, 
we  define 


The  result  of  averaging  the  l  —  2  harmonic  over  the  mean  anomaly  of  the  satellite  yields 

<d)2  =  0 
</>  2  = 


2  2 

3  s3  (e  s2  sin(v)  -  s2  sin(w)  +  4  e  si  cos(b)  +  si  cos(b))  za 


2  1/2 
2  (1  -  •  ) 
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<n}2  = 

2  2 

3  s3  (4  e  si  sin(w)  +  si  sin(y)  -  e  s2  cos(s)  +  s2  cos(b))  za 


2  1/2 

2  (1  -  e  )  sin(i) 


<e}2  = 

2  1/2 

15  e  (1  -  e  )  si  s2  za 


2 


(a.’}  2  =  —  COsI(fl)2  — 

21/22  2 
3  (1  -  e  )  (s2  *  4  si  ♦  i)  za 


2 


(M  —  n0)2  =  —  \  1  —  e2  (u>)2  +  cos  I  {Cl)  ^ 


+ 


2  2  2  2  2  2  2 
(6  e  s2  -  6  s2  -  33  e  si  -  6  si  +  9  e  +  4)  za 


2 


with  short  period  corrections 


nia2  = 

2 

-  a  (6  e  si  s2  sin(2  E)  -  6  si  s2  sin(2  E) 

221/22  2  1/2  2 

-  3  e  (1  -  e  )  s2  cos(2  E)  +  3  (1  -  e  )  s2  cos(2  E) 

21/22  2  2  1/2  3 

-  3  (1  -  e  )  si  cos(2  E)  +  o  (1  -  e  )  cos(2  E)  -  12  c  si  s2  sin(E) 

21/22  2  1/2 

+  12  e  si  s2  sin(E)  +  12  e  (1  ■  e  )  si  cos(E)  -  4  •  (1  -  e  )  cos(E)) 
2  1/2 

/ (2  (1  -  e  )  ) 


za 
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nSL  = 

3  2  1/2 

-  s3  (c  (1  -  e  )  s2  sin(3  E)  sin(w) 

2  1/2  3 

e  (1  -  c  )  s 2  sin(3  E)  sin(w)  +  e  si  cos(3  E)  sin(y) 

2  2  1/2 

c  si  cos(3  E)  sin(y)  -3c  (l  -  c  )  s2  sin(2  E)  sin(w) 

2  1/2  4 

3  (1  -  c  )  s2  sin(2  E)  sin(w)  -  3  e  si  cos(2  E)  sin(w) 

3  2  1/2 

3  si  cos(2  E)  sin(u)  -  3  e  (1  -  e  )  s2  sin(E)  sin(w) 

2  1/2  3 

3  e  (l  -  e  )  s2  sin(E)  sin(w)  +  15  e  si  cos(E)  sin(w) 

2  1/2 

15  e  si  cos(E)  sin(y)  -  e  (1  -  e  )  si  sin(3  E)  cos(s) 

3 

e  s2  cos(3  E)  cos(b)  +  e  s2  cos(3  E)  cos(w) 

221/2  2  1/2 
6  e  (1  -  e  )  si  sin(2  E)  cos(s)  +  3  (1  -  e  )  si  sin(2  E)  cos(w) 
4 

3c  s2  cos(2  E)  cos(v)  -  3  s2  cos(2  E)  cos(b) 

321/2  2  1/2 

12  e  (1  -  c  )  si  sii\(E)  cos(v)  -  33  c  (1  -  e  )  si  sin(E)  cos(v) 
3 

15  c  s2  cos(E)  cos(r)  +  15  e  s2  cos(E)  cos(v))  za/(4  (e  -  1)  (c  +  l)) 


= 

2  1/2  3 

s3  (c  (1  -  c  )  si  sin(3  E)  sin(y)  +  e  s2  cos(3  E)  sin(w) 
2  2  1/2 

c  s2  cos(3  E)  sin(v)  -6c  (l  -  c  )  si  sin(2  E)  sin(w) 

2  1/2  4 

3  (1  -  e  )  si  sin(2  E)  sin(v)  -3c  s2  cos(2  E)  sin(w) 

3  2  1/2 

3  s2  cos(2  E)  sin(B)  +  12  e  (1  -  e  )  si  sin(E)  sin(w) 

2  1/2  3 

33  e  (1  -  e  )  si  sin(E)  sin(y)  ^  15  e  s2  cos(E)  sin(w) 

3  2  1/2 

15  e  s2  cos(E)  sin(s)  +  c  (1  -  e  )  s2  sin(3  E)  cos(y) 

2  1/2  3 

e  (1  -  e  )  s2  sir\(3  E)  cos(w)  +  e  si  cos(3  E)  cos(w) 

2  2  1/2 

e  si  cos(3  E)  cos(b)  -  3  e  (1  -  e  )  s2  sin(2  E)  cos(w) 

2  1/2  4 

3  (1  -  e  )  s2  sin(2  E)  cos(w)  -3c  si  cos(2  E)  cos(b) 

3  2  1/2 

3  si  cos(2  E)  cos(b)  -  3  e  (1  -  e  )  s2  sin(E)  cos(w) 

2  1/2  3 

3  e  (1  -  e  )  s2  sin(E)  cos(b)  +  15  e  si  cos(E)  cos(w) 

15  e  si  cos(E)  cos(w))  za/(4  (c  -  1)  (e  +  1)  sin(i)) 
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nSe2  =  ^n*a2+ 

321/2  2  1/2 

“  (e  (1  -  c  )  si  s2  sin(3  E)  -  2  e  (1  -  e  )  si  s2  sin(3  E) 
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+  3  a  (1  -  e  )  si  s2  sin(2  E)  +  6  (1  -  e  )  si  s2  sin(2  E) 
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-  63  (1  -  a  )  si  sin (2  E)  +  63  a  (1  -  a  )  sin(2  E) 

4  2 

-  72  e  si  s2  cos(2  E)  -  54  a  si  s2  cos(2  E)  +  126  si  s2  cos(2  E) 

321/22  2  1/2  2 

-  9  a  (1  -  a  )  s2  sin(E)  +  9  a  (l  -  a  )  s2  sin(E) 

321/22  2  1/2  2 

+  144  a  (1  -  e  )  si  sin(E)  +  639  a  (1  -  a  )  si  sin(E) 

321/2  2  1/2  4  2 

-  45  a  (l  -  a  )  sin(E)  -  216  e  (1  -  a  )  sin(E)  -  108  a  si  s2  cos  (E) 
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2  2  3 

+  108  ©  si  s2  cos  (E)  +  630  ©  si  s2  cos(E)  -  630  ©  si  s2  cos(E) 
4  2  21/2 

+  54  ©  si  s2  -  54  e  si  s2)  za/(24  (l  -  ©  )  ) 


The  result  of  averaging  the  /  =  3  harmonic  over  the  mean  anomaly  of  the  satellite  yields 


(«}•  =  o 


</>3  = 

2  2  2 
-  15  e  s3  (10  ©  si  s2  sin(v)  -  10  si  s2  sin(w)  -  5  e  s2  cos(w) 

2  2  2  2  2 
+  5  s2  cos(v)  +  20  e  si  cos(e)  +  15  si  cos(w)  -  3  ©  cos(w)  -  4  cos(s)) 
2  1/2 

za2/(l6  (1  -  ©  )  ) 


<6/3  = 

2  2  2  2  2 
15  e  s3  (5  e  s2  sin(y)  -  5  s2  sin(v)  -  20  e  si  sin(v) 
2  2  2 

-  15  si  sin(y)  +  3  ©  sin(v)  +  4  sin(y)  +  10  ©  si  s2  cos(y) 

2  1/2 

-  10  si  s2  cos(v))  2a2/(16  (1  -  ©  )  sin(i)) 


(e}3  = 


2  1/2  22  2  22  2  2 
15  (1  -  ©  )  s2  (5  ©  s2  -  5  s2  -  30  ©  si  -  5  si  +  3  e  +4)  za2 


16 


(u-7 } 3  =  -  COE  I 

2  1/2  22  2  22  2  2 
15  (1  -  ©  )  si  (15  ©  s2  -  5  s2  -  20  ©  si  -  5  si  +  9  ©  +4)  2a2 


16  e 
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M  -  n0)3  =  -\  1  -  e2  [(J;)3  +  cos I(fl)s 


+ 


2  2  2  2  2  2  2 
3  e  sl  (165  e  s2  -  165  s2  -  260  e  sl  -  165  si  +  123  e  +  132)  za2 


16 


with  short  period  corrections 


nSaz  = 

4  3  2  3  3 

-  a  (5  c  s2  sin(3  E)  -  10  e  s2  sin(3  E)  +  5  s2  sin(3  E) 

2  2  2  4 

+  15  c  sl  s2  sin(3  E)  -  15  sl  s2  sin(3  E)  -  3  e  s2  sin(3  E) 

2  221/22 

+  3  c  s2  sin(3  E)  -  15  e  (1  -  e  )  sl  s2  cos(3  E) 

21/22  2  1/2  3 

+  15  (1  •  •  )  sl  s2  cos (3  E)  -  5  (1  -  •  )  sl  cos(3  E) 

221/2  32  2 

+  3  e  (l  -  e  )  sl  cos(3  E)  -  60  e  sl  s2  sin(2  E)  +  60  c  sl  s2  sin(2  E) 

3  321/22 

+  12  e  s2  sin(2  E)  -  12  e  s2  sin(2  E)  +  30  e  (1  -  c  )  sl  s2  cos(2  E) 
21/22  2  1/2  3 

-  30  e  (1  -  «  )  sl  s2  cos(2  E)  +  30  e  (l  -  e  )  sl  cos(2  E) 

321/2  2  1/2 

-  6  c  (1  *  e  )  sl  cos(2  E)  *  12  e  (1  -  e  )  sl  cos(2  E) 

4  3  2  3  3  4  2 

-  15  e  s2  sin(E)  +  30  e  s2  sin(E)  -  15  s2  sin(E)  +  60  e  sl  s2  sin(E) 

2  2  2  4  2 

-  45  e  sl  s2  sin(E)  -  15  sl  s2  sin(E)  -3c  s2  sin(E)  -9c  s2  sin(E) 

2  21/2  2 

+  12  s2  sin(E)  +  15  e  (1  -  e  )  sl  s2  cos(E) 

21/22  22  1/2  3 

-  15  (1  -  e  )  sl  s2  cos(E)  -  60  c  (1  -  e  )  sl  cos(E) 

2  1/2  3  2  2  1/2 

-  15  (1  -  e  )  sl  cos(E)  +  33  e  (1  -  c  )  sl  cos(E) 

2  1/2  2  1/2 
+  12  (1  -  e  )  sl  cos(E))  za2/(4  (1  -  c  )  ) 


nSIz  = 

3 

s3  (30  c  sl  s2  sin(4  E)  sin(w)  -  30  e  sl  s2  sin(4  E)  sin(w) 
321/22  2  1/2  2 

-  15  c  (1  -  e  )  s2  cos(4  E)  sin(w)  +  15  e  (1  -  c  )  s2  cos(4  E)  sin(w) 

21/22  32  1/2 

-  15  e  (1  -  c  )  sl  cos(4  E)  sin(w)  +  3  e  (1  -  e  )  cos(4  E)  sin(w) 

4 

-  80  e  sl  s2  sin(3  E)  sin(w)  +  80  sl  s2  sin(3  E)  sin(w) 

221/22  2  1/2  2 
+  40  c  (1  -  c  )  s2  cos (3  E)  sin(w)  -  40  (1  -  e  )  s2  cos(3  E)  sin(w) 
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221/22  2  1/2  2 
80  e  (1  -  c  )  si  cos(3  E)  sin(w)  +  40  (1  -  e  )  si  cos(3  E)  sin(w) 
221/2  3 

24  e  (1  -  c  )  cos(3  E)  sin(w)  +  240  e  si  s2  sin(2  E)  sin(w) 

3  21/22 

240  e  si  s2  sin(2  E)  sin(s)  +  60  e  (1  -  e  )  s2  cos(2  E)  sin(w) 

2  1/2  2 

60  e  (1  -  c  )  s2  cos(2  E)  sin(w) 

3  21/22 

120  e  (1  -  e  )  si  cos(2  E)  sin(w) 

21/22  32  1/2 

300  e  (1  -  c  )  si  cos(2  E)  sin(s)  +  12  e  (1  -  e  )  cos(2  E)  sin(w) 

2  1/2  4 

72  e  (1  -  c  )  cos(2  E)  sin(s)  +  240  e  si  s2  sin(E)  sin(w) 

2  21/22 

240  si  s2  sin(E)  sin(s)  -  360  e  (1  -  c  )  s2  cos(E)  sin(v) 

21/22  22  1/2  2 

360  (1  -  e  )  s2  cos(E)  sin(w)  +  720  e  (1  -  c  )  si  cos(E)  sin(w) 
21/22  22  1/2 

120  (1  -  c  )  si  cos(E)  sin(w)  -  72  e  (1  -  e  )  cos(E)  sin(w) 

21/2  32 

96  (1  -  e  )  cos(E)  sin(w)  -  15  c  s2  sin(4  E)  cos(w) 

2  2  3 

15  e  s2  sin(4  E)  cos(v)  -  15  e  si  sin(4  E)  cos(v)  +  3  e  sin(4  E)  cos(v) 
21/2  42 

30  «  (1  -  e  )  si  s2  cos(4  E)  cos(p)  +  40  e  s2  sin(3  E)  cos(v) 

2  2  2  2 
40  62  sin(3  E)  cos(v)  +  120  c  si  sin(3  E)  cos(w)  +  40  si  sin(3  E)  cos(w) 
4  2 

8  e  sin(3  E)  cos(v)  -  24  e  sin(3  E)  cos(v) 

2  2  1/2 

160  c  (1  -  c  )  si  s2  cos(3  E)  cos(w) 

21/2  32 

80  (1  -  e  )  si  s2  cos(3  E)  cos(v)  -  120  c  s2  sin(2  E)  cos(v) 

2  3  2 

120  e  s2  6in(2  E)  cos(v)  -  360  e  si  sin(2  E)  cos(w) 

2  3 

480  e  si  sin(2  E)  cos(v)  +  96  e  sin(2  E)  cos(w)  +  72  e  sin(2  E)  cos(w) 

3  2  1/2 

240  e  (1  -  e  )  si  s2  cos(2  E)  cos(w) 

21/2  42 

600  e  (1  *  e  )  si  s2  cos(2  E)  cob(b)  -  120  e  s2  sin(E)  cos(w) 

2  4  2  2  2 

120  s2  sin(E)  cos(s)  +  480  e  si  sin(E)  cos(w)  +  2520  e  si  sin(E)  cos(s) 
2  4  2 

360  si  sin(E)  cos(w)  -  72  c  sin(E)  cos(w)  -  504  e  sin(E)  cos(w) 

2  2  1/2 

96  sin(E)  cos(v)  -  1440  e  (1  -  e  )  si  s2  cos(E)  cos(v) 

2  1/2  2  1/2 
240  (1  “  c  )  si  s2  cos(E)  cos(b))  z&2/(64  (1  -  e  )  ) 


nSQs  — 

3  2  2 

-  s3  (15  e  s2  sin(4  E)  sin (b)  -  15  e  s2  sin(4  E)  sin(v) 
2  3 

15  e  si  sin(4  E)  sin(s)  -  3  e  sin(4  E)  sinCs) 
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21/2  42 

-  30  •  (1  -  •  )  si  s2  cos(4  E)  sin(w)  -  40  e  s2  sin(3  E)  sin(w) 

2  2  2  2 
+  40  s2  sin(3  E)  sin(w)  -  120  ©  si  sin(3  E)  sin(s)  -  40  si  sin(3  E)  sin(w) 
4  2 

+  8  ©  sin(3  E)  sin(w)  +  24  e  sin(3  E)  sin(w) 

2  2  1/2 

+  160  ©  (1  -  e  )  si  s2  cos(3  E)  sin(w) 

21/2  32 

+  80  (1  -  ©  )  si  s2  cos(3  E)  sin(w)  +  120  ©  s2  sin(2  E)  sin(w) 

2  3  2 

-  120  ©  s2  sin(2  E)  sin(w)  +  360  ©  si  sin(2  E)  sin(s) 

2  3 

+  480  ©  si  sin(2  E)  sin(w)  -  96  ©  sin(2  E)  sin(v)  -  72  ©  sin(2  E)  sin(w) 

3  2  1/2 

-  240  ©  (1  -  ©  )  si  s2  cos(2  E)  sin(w) 

21/2  42 

-  600  e  (1  -  e  )  si  s2  cos(2  E)  sin(w)  +  120  e  s2  sin(E)  sin(w) 

2  4  2  2  2 

-  120  s2  sin(E)  sin(w)  -  480  e  si  sin(E)  sin(w)  -  2520  ©  si  sin(E)  sin(w) 

2  4  2 

-  360  si  sin(E)  sin(w)  +  72  ©  sin(E)  sin(w)  +  504  ©  sin(E)  sin(w) 

2  2  1/2 

*►  96  sin(E)  sin(r)  +  1440  ©  (l  -  ©  )  si  s2  cos(E)  sin(w) 

2  1/2  3 

*►  240  (1  -  ©  )  si  s2  cos(E)  sin(w)  +  30  e  si  s2  sin(4  E)  cos(w) 

3  21/22 

-  30  e  si  s2  sin(4  E)  cos(w)  -  15  ©  (1  -  ©  )  s2  cos(4  E)  cos(w) 

21/22  2  1/2  2 
+  15  ©  (1  -  ©  )  s2  cos(4  E)  cos(w)  -  15  ©  (1  -  ©  )  si  cos(4  E)  cos(w) 
321/2  4 

+  3  ©  (1  -  ©  )  cos(4  E)  cos(w)  -  80  ©  si  s2  sin(3  E)  cos(w) 

2  21/22 

+  80  si  s2  sin(3  E)  cos(w)  ♦  40  ©  (1  -  ©  )  s2  cos(3  E)  cos(w) 

21/22  22  1/2  2 

-  40  (1  -  ©  )  s2  cos(3  E)  cos(w)  +  80  ©  (l  -  ©  )  si  cos(3  E)  cos(f) 

21/22  22  1/2 

+  40  (1  -  ©  )  si  cos(3  E)  cos(w)  -  24  ©  (l  -  ©  )  cos(3  E)  cos(w) 

3 

+  240  ©  si  s2  sin(2  E)  cos(w)  -  240  ©  si  s2  sin(2  E)  cos(y) 

321/22  2  1/2  2 

+  60  ©  (1  -  ©  )  s2  cos(2  E)  cos(w)  -  60  ©  (1  -  ©  )  s2  cos(2  E)  cos(w) 

3  21/22 

-  120  e  (1  -  e  )  si  cos(2  E)  cos(w) 

21/22  32  1/2 

-  300  ©  (1  -  ©  )  si  cos(2  E)  cos(w)  +  12  ©  (1  -  ©  )  cos(2  E)  cos(w) 

2  1/2  4 

+  72  ©  (l  -  ©  )  cos(2  E)  cos(v)  +  240  ©  si  s2  sin(E)  cos(w) 

2  21/22 

-  240  si  s2  sin(E)  cos(w)  -  360  ©  (1  -  ©  )  s2  cos(E)  cos(w) 

21/22  22  1/2  2 
+  360  (1  -  e  )  s2  cos(E)  cos(w)  +  720  ©  (1  -  ©  )  si  cos(E)  cos(w) 

21/22  22  1/2 

+  120  (1  -  ©  )  si  cos(E)  cos(w)  -  72  ©  (l  -  ©  )  cos(E)  cos(w) 

2  1/2  2  1/2 

-  96  (1  -  ©  )  cos(E)  cos(w))  za2/(64  (1  -  e  )  sin(i)) 
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n6e3  =  lj£-n6a3~ 

21/233  3 

(1  -  e  )  (15  e  s2  sin(4  E)  -  15  e  62  sin(4  E) 

3  2  2  3 

-  30  e  si  s2  sin(4  E)  +  45  e  si  s2  sin(4  E)  -  3  e  s2  sin(4  E) 

321/22  2  1/2  2 

+  15  e  (1  -  e  )  si  s2  cos(4  E)  -  45  e  (1  -  e  )  el  s2  cos(4  E) 
21/23  3  2  1/2 

+  15  e  (1  -  e  )  si  cos (4  E)  -  3  e  (l  -  e  )  si  cos(4  E) 

4  3  3  4  2 

-  40  e  s2  sin(3  E)  +  40  s2  sin(3  E)  +  80  e  si  s2  sin(3  E) 

2  2  2  4 

-  120  e  si  s2  sin(3  E)  -  120  si  s2  sin(3  E)  +  8  e  s2  sin(3  E) 

2  221/22 
+  24  t  s2  sin(3  E)  +  120  e  (l  -  e  )  si  s2  cos(3  E) 

21/22  22  1/2  3 

+  120  (1  -  e  )  si  s2  cos(3  E)  -  80  e  (l  -  e  )  si  cos(3  E) 
21/23  22  1/2 

-  40  (1  -  c  )  si  cos(3  E)  +  24  e  (1  -  e  )  si  cos(3  E) 

3  3  3  3  2 

+  120  e  s2  sin( 2  E)  -  120  e  s2  sin(2  E)  +  120  e  si  s2  sin(2  E) 

2  3 

+  720  e  si  s2  sin(2  E)  *  96  e  s2  sin(2  E)  -  72  e  s2  sin(2  E) 
321/22  2  1/2  2 

-  300  e  (l  -  e  )  si  s2  cos(2  E)  -  540  e  (1  -  e  )  si  s2  cos(2  E) 

321/23  2  1/2  3 

+  120  e  (1  -  e  )  si  cos(2  E)  +  300  e  (1  -  e  )  si  cos(2  E) 
321/2  2  1/2 

-  12  e  (l  -  e  )  si  cos(2  E)  -  72  e  (1  -  e  )  si  cos(2  E) 

4  3  3  4  2 

+  120  e  s2  sin(E)  -  120  s2  sin(E)  -  720  e  si  s2  sin(E) 

2  2  2  4 

-  2520  e  si  s2  sin(E)  -  120  si  s2  sin(E)  +  72  e  s2  sin(E) 

2  221/22 
+  504  e  s2  sin(E)  +  96  s2  sin(E)  +  1800  c  (l  -  e  )  si  s2  cos(E) 
21/22  22  1/2  3 

-  120  (1  -  e  )  si  s2  cos(E)  -  720  e  (l  -  e  )  si  cos(E) 

21/23  2  2  1/2 

-  120  (1  *  e  )  si  cos(E)  +  72  e  (1  -  e  )  si  cos(E) 

2  1/2 

+  96  (1  -  c  )  si  cos(E))  za2/(64  e) 


nSuJs  —  —  cos(I)n6Q3  — 

4  2  2  2  2 

-  (15  e  si  s2  sin(4  E)  -  60  e  si  s2  sin(4  E)  +  45  si  s2  sin(4  E) 

2  3  3  4  2 

+  15  e  si  sin(4  E)  -  15  si  sin(4  E)  -  3  e  si  sin(4  E)  +  3  e  si  sin(4  E) 
221/23  2  1/2  3 

+  15  e  (1  -  e  )  s2  cos(4  E)  -  15  (1  -  e  )  s2  cos(4  E) 
221/22  2  1/2  2 

-  30  e  (l  *  e  )  si  s2  cos(4  E)  +  45  (1  -  e  )  si  s2  cos(4  E) 

221/2  52  32 

-  3  e  (l  -  e  )  s2  cos(4  E)  +  40  e  si  s2  sin(3  E)  +  80  e  si  s2  sin(3  E) 
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2  3  3  3 

-  120  e  si  s2  sin(3  E)  -  40  e  si  sin(3  E)  +  40  e  si  sin(3  E) 

5  3  321/23 

-8c  si  sin(3  E)  +  8  a  si  sin(3  E)  -  40  c  (1  -  e  )  s2  cos(3  E) 

21/23  2  1/2  2 

+  40  e  (1  -  e  )  s2  cos(3  E)  -  120  e  (1  -  e  )  si  s2  cos(3  E) 
321/2  42 

+  24  e  (l  -  e  )  s2  cos(3  E)  -  480  e  si  s2  sin(2  E) 

2  2  2  4  3 

+  720  e  si  s2  sin(2  E)  -  240  si  s2  sin(2  E)  -  120  e  si  sin(2  E) 

3  4  2 

+  120  si  sin(2  E)  *  168  e  si  sin(2  E)  -  144  c  si  sin(2  E)  -  24  si  sin(2  E) 
221/23  2  1/2  3 

-  60  c  (1  -  e  )  s2  cos(2  E)  +  60  (1  -  e  )  s2  cos(2  E) 

4  21/22  2  21/22 

+  240  e  (1  -  e  )  si  s2  cos(2  E)  +  480  e  (1  -  a  )  si  s2  cos(2  E) 

21/22  42  1/2 

-  300  (1  -  e  )  si  s2  cos(2  E)  -  48  e  (1  -  e  )  s2  cos(2  E) 

221/2  2  1/2 

-  60  e  (l  -  e  )  *  s2  cos(2  E)  +  24  (1  -  e  )  s2  cos(2  E) 

5  2  3  2  2 

-  120  e  si  s2  sin(E)  -  240  e  si  s2  sin(E)  +  360  a  si  s2  sin(E) 

5  3  3  3  3 

+  480  e  si  sin(E)  +  1080  e  si  sin(E)  -  1560  e  si  sin(E) 

5  3 

-  264  e  si  sin(E)  -  600  e  si  sin(E)  +  864  a  si  sin(E) 

321/23  2  1/2  3 

+  360  e  (1  -  e  )  s2  cos(E)  -  360  a  (1  -  a  )  s2  cos(E) 
321/22  2  1/2  2 

-  2400  a  (1  -  e  )  si  s2  cos(E)  +  1560  e  (1  -  a  )  si  s2  cos(E) 

321/2  2  1/2 

+  264  a  (1  -  e  )  s2  cos(E)  -  96  a  (1  -  a  )  s2  cos(E))  za2 
2  1/2 

/ (64  a  (1  -  c  )  ) 


n6M3  =  — v  1  —  c2  jiSu.'3  4-  cos (I)nSQ3]  4- 

321/2  3  2  1/2  3 

3  (-  15  a  (1  -  a  )  cos(4  E)  s2  +  15  a  (l  -  a  )  cos(4  E)  s2 
221/2  3  2  1/2  3 

+  40  a  (1  -  e  )  cos(3  E)  s2  -  40  (1  -  a  )  cos(3  E)  s2 

321/2  3  2  1/2  3 

+  60  e  (1  -  a  )  cos(2  E)  s2  -  60  a  (1  -  a  )  cos(2  E)  s2 
221/2  3  2  1/2  3 

-  360  e  (1  -  a  )  cos(E)  s2  +  360  (1  -  a  )  cos(E)  s2 

3  2  2  4  2 

+  45  a  sin(4  E)  si  s2  -  45  a  sin(4  E)  si  s2  -  120  a  sin(3  E)  si  s2 
2  3  2  2 

+  120  sin(3  E)  si  s2  +  360  a  sin(2  E)  si  s2  -  360  a  sin(2  E)  si  s2 
4  2  2  2  2 

+  120  a  sin(E)  si  s2  +  240  a  sin(E)  si  s2  -  360  sin(E)  si  s2 

21/2  2  22  1/2  2 

-  45  a  (1  -  a  )  cos(4  E)  si  s2  +  240  a  (1  -  a  )  cos(3  E)  si  s2 

21/2  2  32  1/2  2 

+  120  (1  -  a  )  cos(3  E)  si  s2  -  360  a  (1  -  a  )  cos(2  E)  si  s2 
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21/2  2  22  1/2  2 

-  900  e  (1  -  e  )  cos(2  E)  si  s2  +  2160  e  (1  -  t  )  cos(E)  si  s2 

21/2  2  32  1/2 

+  360  (1  -  e  )  cos(E)  si  s2  +  9  e  (l  -  e  )  cos(4  E)  s2 
221/2  32  1/2 

-  72  c  (1  -  c  )  cos(3  E)  s2  +  36  e  (1  -  e  )  cos(2  E)  s2 

21/2  22  1/2 
+  216  e  (1  -  e  )  cos(2  E)  s2  -  216  c  (1  -  «  )  cos(E)  s2 

21/2  323 

-  288  (1  -  e  )  cos(E)  s2  +  15  e  sin(4  E)  si  -  120  e  sin(3  E)  si 

3  3  3  3 

-  40  sin(3  E)  si  +  360  e  sin(2  E)  si  +  480  e  sin(2  E)  si 

4  3  2  3  3  3 

-  320  c  sin(E)  si  -  2280  c  sin(E)  si  -  360  sin(E)  si  -  9  •  sin(4  E)  si 

4  2  3 

+  24  e  sin(3  E)  si  +  72  c  sin(3  E)  si  -  288  e  sin(2  E)  si 

4  2 

-  216  e  sin(2  E)  si  +  168  e  sin(E)  si  +  1320  e  sin(E)  si  +  288  sin(E)  si) 
za2/64 


In  order  to  express  the  long  period  averaging  of  the  above  over  the  mean  anomaly  of  the  Moon, 
it  suffices  to  average  various  products  of  51,52  and  53  with  za  for  the  /  =  2  harmonic  and  with 
z a2  for  the  /  =  3  harmonic.  For  example,  if  one  would  like  to  propagate  the  eccentricity  as  it  is 
influenced  by  the  /  =  2  harmonic,  one  uses  the  expressions  for  { e)2  and  be 2  above  along  with  the 
expressions  for  (zaS\S2)  and  6(2a5l52)  below.  After  these  are  evaluated,  the  doubly  averaged 
secular  rate  is  computed 


(<e)>2 


(zaSlS2) 


and  e  is  propagated  by  the  formula 


1  ,  1  _ 

e(t)  =  c0  -r  ((e))2(i  -  to) - ^ - <F(ca5l52)  -  Sc2. 

where  ({e))2  is  a  function  of  the  initial  conditions,  the  long  period  correction  is  a  function  of  the 
updated  mean  elements,  and  be 2  is  a  function  of  the  updated  mean  elements  with  long  period 
corrections. 

The  presentation  of  the  expressions  below  requires  the  definition  of  two  quantities:  for  the 
/  =  2  harmonic  we  define 

-  f2?) 
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and  for  the  /  =  3  harmonic  we  define 


zb2  =  zb(  —  )  . 

\ak/ 

The  result  of  averaging  the  pertinent  products  for  the  rates  from  l  =  2  harmonic  over  the  period 
of  the  Moon  is 

(za)  =  nzb 

(zaSl7)  = 

2  2  2 
(sm  (apk)  cos  (jk)  +  cos  (&pk))  n  2b 


(,a52!)  = 


2  2  2 
(cos  (apk)  cos  (jk)  +  sin  (apk))  n  2b 


(zaSlS2)  = 


(zaSlS3)  = 


[zaS2SZ)  = 


cos(apk)  sin(apk)  sin  (jk)  n  2b 
2 


sin(apk)  cos(jk)  sin(jk)  n  2b 
2 


cos(apk)  cos(jk)  sin(jk)  n  zb 
2 


with  long  period  corrections 
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6(:a)  = 


ek  sin(fk)  n  zb 


nk 


6(cg512)  = 

2  2 

n  (cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(wk)  sin(wk) 

2 

-  ek  cos(3  fk)  cos  (jk)  cos(ik)  sin(wk) 

2  2 

+  3  cos  (apk)  cos(2  fk)  cos  (jk)  cos(wk)  sin(wk) 

2 

-  3  cos(2  fk)  cos  (jk)  cos(wk)  sin(vk) 

2  2 

*+■  3  cos  (apk)  ek  cos(fk)  cos  (jk)  cos(wk)  sin(wk) 

2 

-  3  ek  cos(fk)  cos  (jk)  cos(wk)  sin(wk) 

+  4  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk)  sin(vk) 

+  12  cos(apk)  sin(apk)  cos(fk)  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

+  8  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

2 

+  cos  (apk)  ek  cos(3  fk)  cos(rk)  sin(wk) 

2 

+  3  cos  (apk)  cos(2  fk)  cos(wk)  sin(wk) 

2 

+  3  cos  (apk)  ek  cos(fk)  cos(yk)  sin(wk) 

2  2  2 

2  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2 

-  2  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  2 

+  6  cos  (apk)  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

2  2  2  2 

-  6  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  +  4  cos  (apk)  ek  sin(fk)  cos  (jk) 

2  2  2 

cos  (wk)  -  4  ek  sin(fk)  cos  (jk)  cos  (wk) 

2 

-  2  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 

2 

-  6  cos(apk)  sin(apk)  cos(2  fk)  cos(jk)  cos  (wk) 

2 

-  6  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

2  2 
+  2  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (wk) 

2  2  2  2 
+  6  cos  (apk)  cos(fk)  sin(fk)  cos  (wk)  +  4  cos  (apk)  ek  sin(fk)  cos  (wk) 

2  2  2 

-  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  +  ek  sin(fk)  cos(2  fk)  cos  (jk) 

2  2  2 

-  3  cos  (apk)  cos(fk)  sin(fk)  cos  (jk)  +  3  cos(fk)  sin(fk)  cos  (jk) 
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2  2  2 

-  5  cos  (apk)  ek  sin(fk)  cos  (jk)  +  5  ek  sin(fk)  cos  (jk) 

+  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk) 

+  3  cos(apk)  sin(apk)  cos(2  fk)  cos(jk) 

2 

+  3  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  -  cos  (apk)  ek  sin(fk)  cos(2  fk) 
2  2 

-  3  cos  (apk)  cos(fk)  sin(fk)  +  cos  (apk)  ek  sin(fk))  zb/(6  nk) 


6(zaS  22)  = 

2  2 

-  n  (cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(wk)  sin(yk) 

2  2 

+  3  cos  (apk)  cos(2  fk)  cos  (jk)  cos(yk)  sin(yk) 

2  2 

+  3  cos  (apk)  ek  cos(fk)  cos  (jk)  cos(wk)  sin(yk) 

+  4  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(uk)  sin(yk) 

+  12  cos(apk)  sin(apk)  cos(fk)  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

+  8  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

2 

+  cos  (apk)  ek  cos(3  fk)  cos(yk)  sin(yk)  -  ek  cos(3  fk)  cos(yk)  sin(yk) 

2 

+  3  cos  (apk)  cos (2  fk)  cos(yk)  sin(yk)  -  3  cos(2  fk)  cos(yk)  sin(yk) 

2 

+  3  cos  (apk)  ek  cos(fk)  cos(yk)  sin(yk)  -  3  ek  cos(fk)  cos(wk)  sin(yk) 

2  2  2 
+  2  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (yk) 

2  2  2 
+  6  cos  (apk)  cos(fk)  sm(fk)  cos  (jk)  cos  (yk) 

2  2  2 
+  4  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (yk) 

2 

-  2  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (yk) 

2 

*  6  cos(apk)  sin(apk)  cos(2  fk)  cos(jk)  cos  (yk) 

2 

-  6  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (yk) 

2  2  2 
+  2  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (yk)  -  2  ek  sin(fk)  cos(2  fk)  cos  (yk) 
2  2  2 
+  6  cos  (apk)  cos(fk)  sin(fk)  cos  (yk)  -  6  cos(fk)  sin(fk)  cos  (yk) 

2  2  2 
+  4  cos  (apk)  ek  sm(fk)  cos  (yk)  -  4  ek  sm(fk)  cos  (yk) 

2  2 

-  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk) 

2  2  2  2 

-  3  cos  (apk)  cos(fk)  sin(fk)  cos  (jk)  -  5  cos  (apk)  ek  sin(fk)  cos  (jk) 

+  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk) 

+  3  cos(apk)  sin(apk)  cos(2  fk)  cos(jk) 
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2 

+  3  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  -  cos  (apk)  ek  sin(fk)  cos(2  fk) 

2 

+  ek  sin(fk)  cos(2  fk)  -  3  cos  (apk)  cos(fk)  sin(fk)  +  3  cos(fk)  sin(fk) 

2 

+  cos  (apk)  ek  sin(fk)  -  ek  sin(fk))  zb/(6  nk) 


<5(;a5l52)  = 

2 

-  n  (2  cos(apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos(yk)  sin(yk) 

2 

+  6  cos(apk)  sin(apk)  cos(2  fk)  cos  (jk)  cos(wk)  sin(yk) 

2 

+  6  cos(apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos(wk)  sin(wk) 

2 

-  8  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk)  sin(yk) 

+  4  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(vk)  sin(wk) 

2 

-  24  cos  (apk)  cos(fk)  sin(fk)  cos(jk)  cos(vk)  sin(wk) 

12  cos(fk)  sin(fk)  cos(jk)  cos(yk)  sin(yk) 

2 

-  16  cos  (apk)  ek  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

+  8  ek  sin(fk)  cos(jk)  cos(wk)  sin(wk) 

*♦*  2  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(wk)  sin(sk) 

^  6  cos(apk)  sin(apk)  cos(2  fk)  cos(yk)  sin(yk) 

+  6  cos(apk)  sin(apk)  ek  cos(fk)  cos(yk)  sin(yk) 

2  2 

+  4  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (yk) 

2  2 

+  12  cos(apk)  sin(apk)  cos(fk)  sin(fk)  cos  (jk)  cos  (yk) 

2  2 

+  8  cos(apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (yk) 

2  2  2 
+  4  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos  (yk)  -  2  ek  cos(3  fk)  cos(jk)  cos  (wk) 
2  2  2 
+  12  cos  (apk)  cos(2  fk)  cos(jk)  cos  (yk)  -  6  cos(2  fk)  cos(jk)  cos  (yk) 

2  2  2 
+  12  cos  (apk)  ek  cos(fk)  cos(jk)  cos  (wk)  -  6  ek  cos(fk)  cos(jk)  cos  (yk) 

2 

+  4  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (yk) 

2 

+  12  cos(apk)  sin(apk)  cos(fk)  sin(fk)  cos  (yk) 

2 

+  8  cos(apk)  sin(apk)  ek  sin(fk)  cos  (yk) 

2 

-  2  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk) 

2 

-  6  cos(apk)  sin(apk)  cos(fk)  sin(fk)  cos  (jk) 

2  2 

-  10  cos(apk)  sin(apk)  ek  sin(fk)  cos  (jk)  -  2  cos  (apk)  ek  cos(3  fk)  cos(jk) 
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2 

+  ek  cos(3  fk)  cos(jk)  -  6  cos  (apk)  cos(2  fk)  cos(jk)  +  3  cos(2  fk)  cos(jk) 
2 

-  6  cos  (apk)  ek  cos(fk)  cos(jk)  +  3  ek  cos(fk)  cos(jk) 

•  2  cos (apk)  sin(apk)  ek  sin(fk)  cos(2  fk) 

-  6  cos(apk)  sin(apk)  cos(fk)  sin(fk)  +  2  cos(apk)  sin(apk)  ek  sin(fk))  zt 
/  ( 12  nk) 


*(.aSlS3)  = 

sin(jk)  n  (2  sin(apk)  ek  cos(3  fk)  cos(jk)  cos(wk)  sin(wk) 

+  6  sin(apk)  cos(2  fk)  cos(jk)  cos(wk)  sin(wk) 

+  6  sin(apk)  ek  cos(fk)  cos(jk)  cos(wk)  sin(wk) 

-  4  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(wk)  sin(wk) 

-  12  cos(apk)  cos(fk)  sin(fk)  cos(wk)  sin(wk) 

-  8  cos(apk)  ek  sin(fk)  cos(wk)  sin(wk) 

2 

+  4  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 

2 

+  12  sin(apk)  cos(fk)  sin(fk)  cos(jk)  cos  (wk) 

2  2 

+  8  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk)  +  2  cos(apk)  ek  cos(3  fk)  cos  (wk) 

2  2 
+  6  cos(apk)  cos(2  fk)  cos  (wk)  +  6  cos(apk)  ek  cos(fk)  cos  (wk) 

-  2  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  -  6  sin(apk)  cos(fk)  sin(fk)  cos(jk) 

-  10  sin(apk)  ek  sin(fk)  cos(jk)  -  cos(apk)  ek  cos(3  fk) 

-  3  cos(apk)  cos(2  fk)  -  3  cos(apk)  ek  cos(fk))  zb/(12  nk) 


6(raS2S3)  = 

sin(jk)  n  (2  cos(apk)  ek  cos(3  fk)  cos(jk)  cos(wk)  sin(wk) 
+  6  cos(apk)  cos(2  fk)  cos(jk)  cos(wk)  sin(wk) 

+  6  cos(apk)  ek  cos(fk)  cos(jk)  cos(wk)  sin(wk) 

+  4  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(wk)  sin(wk) 

+  12  sin(apk)  cos(fk)  sin(fk)  cos(wk)  sin(wk) 

+  8  sin(apk)  ek  sin(fk)  cos(wk)  sin(wk) 

2 

+  4  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 
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2 


+  12  cos(apk)  cos(fk)  sin(fk)  cos(jk)  cos  (wk) 

2  2 

+  8  cos(apk)  ek  sin(fk)  cos(jk)  cos  (sk)  -  2  sin(apk)  ek  cos(3  fk)  cos  (wk) 

2  2 

-  6  sin(apk)  cos(2  fk)  cos  (wk)  -  6  sin(apk)  ek  cos(fk)  cos  (wk) 

*  2  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  -  6  cos(apk)  cos(fk)  sin(fk)  cos(jk) 

-  10  cos(apk)  ek  sin(fk)  cos(jk)  +  sin(apk)  ek  cos(3  fk) 

+  3  sin(apk)  cos(2  fk)  ♦  3  sin(apk)  ek  cos(fk))  zb/(12  nk) 


The  result  of  averaging  the  pertinent  products  for  the  rates  from  the  l  =  3  harmonic  over  the  period 
of  the  Moon  is 


-  ek  n  (sin(apk)  cos(jk)  sin(wk)  +  cos(apk)  cos(wk))  zb2 


-  ek  n  (cos(apk)  cos(jk)  sin(vk)  -  sin(apk)  cos(wk))  zb2 


ek  sin(jk)  n  sin(wk)  zb2 


(ra2Sl3)  = 


2  2  2  2 
3  ek  (cos  (apk)  cos  (jk)  -  cos  (jk)  -  cos  (apk))  n 

(sin(apk)  cos(jk)  sin(wk)  +  cos(apk)  cos(wk))  zb2/4 


2 


2 


2 


-  3  ek  (cos  (apk)  cos  (jk)  -  cos  (apk)  +  1)  n 

(cos(apk)  cos(jk)  sin(wk)  -  sin(apk)  cos(wk))  zb2/4 
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(;a2SlS22)  = 

2  3 

-  ek  n  (3  cos  (apk)  sin(apk)  cos  (jk)  sin(wk) 

2 

-  3  cos  (apk)  sin (apk)  cos(jk)  sin(wk)  +  sin(apk)  cos(jk)  sin(wk) 

3  2  2 

+  3  cos  (apk)  cos  (jk)  cos(yk)  -  2  cos(apk)  cos  (jk)  cos(wk) 

3 

-  3  cos  (apk)  cos(wk)  +  3  cos(apk)  cos(wk))  zb2/4 


(za2525l2)  = 

3  3  3 

ek  n  (3  cos  (apk)  cos  (jk)  sin(vk)  -  3  cos(apk)  cos  (jk)  sin(wk) 
3 

*  3  cos  (apk)  cos(jk)  sin(yk)  +  2  cos(apk)  cos(jk)  sin(wk) 

2  2  2 
-  3  cos  (apk)  sin(apk)  cos  (jk)  cos(wk)  +  sin(apk)  cos  (jk)  cos(wk) 

2 

+  3  cos  (apk)  sin(apk)  cos(vk))  zb2/4 


(ca25l5253)  = 

2 

ek  sin(jk)  n  (3  cos(apk)  sin(apk)  cos  (jk)  sin(wk) 

2 

-  cos(apk)  sin(apk)  sin(wk)  +  2  cos  (apk)  cos(jk)  cos(wk)  -  cos(jk)  cos(wk)) 
zb2/4 


(ra2  5  3  522)  = 

2  2  2 

ek  sin(jk)  n  (3  cos  (apk)  cos  (jk)  sin(wk)  -  cos  (apk)  sin(wk)  +  sin(wk) 
-  2  cos(apk)  sin(apk)  cos(jk)  cos(ffk))  zb2/4 


(ra2535l2)  = 

2  2  2 

-  ek  sin(jk)  n  (3  cos  (apk)  cos  (jk)  sin(wk)  -  3  cos  (jk)  sin(wk) 

2 

-  cos  (apk)  sin(wk)  -  2  cos(apk)  sin(apk)  cos(jk)  cos(wk))  zb2/4 


with  long  period  corrections 
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*(;fl2Sl)  = 

2 

-  n  (sin(apk)  ek  sin(3  fk)  cos(jk)  sin(wk) 

+  6  sin(apk)  ek  sin(2  fk)  cos(jk)  sin(wk) 

2 

+  9  sin(apk)  ek  sin(fk)  cos(jk)  sin(wk)  +  12  sin(apk) 
2 

+  cos(apk)  ck  cos(3  fk)  sin(wk)  +  6  cos(apk)  ek  cos(2 
2 

*►  3  cos(apk)  ek  cos(fk)  sin(sk)  +  12  cos(apk)  cos(fk) 
2 

-  sin(apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

-  6  sin(apk)  ek  cos(2  fk)  cos(jk)  cos(vk) 

2 

-  3  sin(apk)  ek  cos(fk)  cos(jk)  cos(wk)  -  12  sin(apk) 

2 

+  cos(apk)  ek  sin(3  fk)  cos(ck)  +  6  cos(apk)  ek  sin(2 
2 

+  9  cos(apk)  ek  sin(fk)  cos(wk)  +  12  cos(apk)  sin(fk) 


£(ra252)  = 

2 

-  n  (cos(apk)  ek  sin(3  fk)  cos(jk)  sin(wk) 

6  cos(apk)  ek  sin(2  fk)  cos(jk)  sin(wk) 

2 

♦  9  cos(apk)  ek  sin(fk)  cos(jk)  sin(uk)  +  12  cos(apk) 

2 

-  sin(apk)  ek  cos(3  fk)  sin(wk)  -  6  sin(apk)  ek  cos(2 

2 

-  3  sin(apk)  ek  cos(fk)  sin(sk)  -  12  sin(apk)  cos(fk) 

2 

-  cos(apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

-  6  cos(apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

2 

-  3  cos(apk)  ek  cos(fk)  cos(jk)  cos(wk)  -  12  cos(apk) 

2 

-  sin(apk)  ek  sin(3  fk)  cos(sk)  -  6  sin(apk)  ek  sin(2 

2 

*  9  sin(apk)  ek  sin(fk)  cos(wk)  -  12  sin(apk)  sin(fk) 


Bin(fk)  cos(jk)  sin(wk) 
fk)  sin(wk) 
sin (wk) 


cos(fk)  cos(jk)  cos(wk) 

fk)  cos(wk) 

cos  (trk) )  zb2/ (12  nk) 


sin(fk)  cos(jk)  sin(sk) 

fk)  sin(wk) 

sin(wk) 

cos(fk)  cos(jk)  cos(wk) 
fk)  cos(wk) 
cos(wk))  zb2/(l2  nk) 


£(ra253)  = 

2 

sin(jk)  n  (ek  sin(3  fk)  sin(¥k)  +  6  ek  sin(2  fk)  sin(wk) 

2  2 
+  9  ek  sin(fk)  sin(wk)  +  12  sin(fk)  sin(wk)  -  ek  cos(3  fk)  cos(wk) 

2 

-  6  ek  cos(2  fk)  cos(wk)  -  3  ek  cos(fk)  cos(wk)  -  12  cos(fk)  cos(wk))  zb2 
/ ( 12  nk) 
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«^(ra25l3)  = 

2  2  3  2 

-  n  (24  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  3  2 

sin(wk)  -  24  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  120  cos  (apk)  sin  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2 

-  120  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  3  2 

+  104  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 
2  3  2 

-  104  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  160  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(vk) 

3  2 

-  160  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  360  cos  (apk)  sin(apk)  ek  cos(fk)  sm(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2 

-  360  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  3  2 

+  112  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

-  112  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  80  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2 

-  80  sin(apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  36  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  36  cos(apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  180  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  180  cos(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  120  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  120  cos(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  240  cos  (apk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  240  cos(apk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  360  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  360  cos(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  180  cos  (apk)  ek  cos(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  180  cos(apk)  ek  cos(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  72  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  360  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 
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2  2  2 
+  312  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 
2  2 
+  480  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  1080  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  336  cos  (apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  240  cos  (apk)  sin(apk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2  2 

+  12  cos  (apk)  ek  cos(5  fk)  cos  (wk)  sin(wk) 

3  2 

+  60  cos  (apk)  ek  cos(4  fk)  cos  (wk)  sin(wk) 

3  2  2 

+  40  cos  (apk)  ek  cos(3  fk)  cos  (wk)  sin(wk) 

3  2 

+  80  cos  (apk)  cos(3  fk)  cos  (wk)  sin(wk) 

3  2 

+  120  cos  (apk)  ek  cos(2  fk)  cos  (wk)  sin(wk) 

3  2  2 

+  60  cos  (apk)  ek  cos(fk)  cos  (wk)  sin(wk) 

2  2  3 

-  6  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 

2  3 

+  6  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  ( j k )  sin(wk) 

2  3 

-  30  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

3 

+  30  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2  3 

-  56  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  3 

+  56  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  3 

-  40  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

3 

+  40  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  3 

-  270  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

3 

*  270  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  3 

-  178  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2  3 

+  178  sin(apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2  3 

-  200  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  sin(wk) 

3 

+  200  sin(apk)  sin(fk)  cos  (jk)  sin(wk) 

3  2  2 

-  9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

2  2 

+  9  cos(apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

3  2 

-  45  cos  (apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

2 

+  45  cos(apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 
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3  2  2 

-  45  cos  (apk)  ek  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 

+  45  cos(apk)  ek  cos(3  fk)  cos  (jk)  sin(wk) 

3  2 

-  60  cos  (apk)  cos(3  fk)  cos  (jk)  sin(wk) 

2 

+  60  cos(apk)  cos (3  fk)  cos  (jk)  sin(wk) 

3  2 

-  180  cos  (apk)  ck  cos(2  fk)  cos  (jk)  sin(wk) 

2 

+  180  cos(apk)  ek  cos(2  fk)  cos  (jk)  sin(wk) 

3  2  2 

-  90  cos  (apk)  ek  cos(fk)  cos  (jk)  sin(wk) 

2  2 

+  90  cos(apk)  ek  cos(fk)  cos  (jk)  sin(wk) 

3  2 

-  180  cos  (apk)  cos(fk)  cos  (jk)  sin(wk) 

2 

+  180  cos(apk)  cos(fk)  cos  (jk)  sin(wk) 

2  2 

-  18  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(wk) 

2 

-  90  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

2  2 

-  48  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

2 

-  120  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

2 

-  90  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(wk) 

2  2 

+  66  cos  (apk)  sin(apk)  ek  sin(fk)  cos(jk)  sin(wk) 

2 

+  120  cos  (apk)  sin(apk)  sin(fk)  cos(jk)  sin(wk) 

3  2  3 

-  3  cos  (apk)  ek  cos(5  fk)  sin(wk)  -  15  cos  (apk)  ek  cos(4  fk)  sin(wk) 

3  2  3 

+  5  cos  (apk)  ek  cos(3  fk)  sin(wk)  -  20  cos  (apk)  cos(3  fk)  sin(wk) 

3  3  2 

+  60  cos  (apk)  ek  cos(2  fk)  sin(wk)  +  30  cos  (apk)  ek  cos(fk)  sin(wk) 

3  2  2  3 

+  180  cos  (apk)  cos(fk)  sin(wk)  -  12  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk) 
3  2  3  3 

cos  (wk)  +  12  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk) 

2  3  3 

-  60  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk) 

3  3 

+  60  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3  3 

-  40  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

+  40  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

-  80  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3 

+  80  sin(apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

-  120  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 
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3  3 

+  120  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3  3 

-  60  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  3  3 

+  60  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

3  2  2  3 

+  72  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3 

,  -  72  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

3  2  3 

+  360  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3 

-  360  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  2  2  3 

+  312  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  312  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3 

+  480  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  3 

-  480  cos(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3 

+  1080  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

2  3 

-  1080  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

3  2  2  3 

+  336  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  336  cos(apk)  ek  sin(fk)  cos  (jk)  cos  (wk) 

3  2  3 

+  240  cos  (apk)  sin(fk)  cos  (jk)  cos  (wk) 

2  3 

-  240  cos(apk)  sin(fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  36  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk) 

2  3 

-  180  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk) 

2  2  3 

-  120  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

-  240  cos  (apk)  sin(apk)  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

-  360  cos  (apk)  sm(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk) 

2  2  3 

-  180  cos  (apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

3  2  3 

+  24  cos  (apk)  ek  sin(fk)  cos(4  fk)  co6  (wk) 

'  3  3 

+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (wk) 

3  2  3 

+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (wk) 

•  3  3 

+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (wk) 

3  3 

+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (wk) 

3  2  3  3  3 

+  112  cos  (apk)  ek  sin(fk)  cos  (wk)  +  80  cos  (apk)  sin(fk)  cos  (wk) 
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2  2  3 

+  9  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos(sk) 

2  3 

-  9  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos(sk) 

2  3 

+  45  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

3 

-  45  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

2  2  3 

+  45  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 
2  3 

-  45  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

2  3 

+  60  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos(wk) 

3 

-  60  sin(apk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  3 

+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

3 

*  180  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

2  2  3 

+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

2  3 

-  90  sin(apk)  ek  cos(fk)  cos  (jk)  cos(yk) 

2  3 

+  180  cos  (apk)  sin(apk)  cos(fk)  cos  (jk)  cos(wk) 

3 

-  180  sin(apk)  cos(fk)  cos  (jk)  cos(wk) 

3  2  2 

-  54  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(¥k) 

2  2 

+  54  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

3  2 

-  270  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(tk) 

2 

+  270  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(wk) 

3  2  2 

-  264  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 

+  264  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(yk) 

3  2 

-  360  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2 

+  360  cos(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

3  2 

-  990  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(wk) 

2 

+  990  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(wk) 

3  2  2 

-  402  cos  (apk)  ek  sin(fk)  cos  (jk)  cos(wk) 

2  2 

+  402  cos(apk)  ek  sin(fk)  cos  (jk)  cos(wk) 

3  2 

-  360  cos  (apk)  sin(fk)  cos  (jk)  cos(wk) 

2 

+  360  cos(apk)  sin(fk)  cos  (jk)  cos(wk) 

2  2 

+  27  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 
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2 

+  135  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

2  2 

+  75  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  180  cos  (apk)  sin(apk)  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

2  2 

+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos(wk) 

2 

-  180  cos  (apk)  sin(apk)  cos(fk)  cos(jk)  cos(wk) 

3  2 

-  18  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(wk) 

3 

-  90  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(wk) 

3  2 

-  48  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(wk) 

3 

-  120  cos  (apk)  sin(fk)  cos(2  fk)  cos(wk) 

3  3  2 

-  90  cos  (apk)  ek  cos(fk)  sin(fk)  cos(wk)  +  66  cos  (apk)  ek  sin(fk)  cos(wk) 

3 

+  120  cos  (apk)  sin(fk)  cos(vk))  zb2/(240  nk) 


6(za2S23)  = 


3  2  3  2 

n  (24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 
3  3  2 

+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  3  2 

+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  («k)  sin(sk) 

3  3  2 

+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  («k)  sm(sk) 

3  3  2 

+  360  cos  (apk)  ek  cos(fk)  siri(fk)  cos  (jk)  cos  (wk)  sin(sk) 

3  2  3  2 

+  112  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  3  2 

+  80  cos  (apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(sk) 

2  2  2  2 

-  36  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  180  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 

-  120  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  240  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  360  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 

-  180  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  72  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 
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2  2 

-  72  cos(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2 

+  360  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  360  cos(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2  2 

+  312  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 
2  2 

-  312  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2 

+  480  cos  (apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  480  cos(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2 

+  1080  cos  (apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  1080  cos(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2  2 

+  336  cos  (apk)  ek  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  336  cos(apk)  ek  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2 

+  240  cos  (apk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  240  cos(apk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 

-  12  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (wk)  sin(wk) 

2  2 

+  12  sin(apk)  ek  cos(5  fk)  cos  (wk)  sin(wk) 

2  2 

-  60  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (wk)  sin(wk) 

2 

+  60  sin(apk)  ek  cos(4  fk)  cos  (wk)  sin(wk) 

2  2  2 

-  40  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (wk)  sin(wk) 

2  2 

+  40  sin(apk)  ek  cos(3  fk)  cos  (wk)  sin(wk) 

2  2 

-  80  cos  (apk)  sin(apk)  cos(3  fk)  cos  (wk)  sin(wk) 

2 

+  80  sin(apk)  cos(3  fk)  cos  (wk)  sin(wk) 

2  2 

-  120  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (wk)  sin(wk) 

2 

+  120  sin(apk)  ek  cos(2  fk)  cos  (wk)  sin(wk) 

2  2  2 

-  60  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (wk)  sin(wk) 

2  2 

+  60  sirv(apk)  ek  cos(fk)  cos  (wk)  sin(wk) 

3  2  3 

-  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 

3  3 

-  30  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

3  2  3 

-  56  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

3  3 

-  40  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 
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3  3 

-  270  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(sk) 

3  2  3 

-  178  cos  (apk)  ek  sin(fk)  cos  (jk)  sin(sk) 

3  3 

-  200  cos  (apk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  2 
+  9  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

2  2 
+  45  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

2  2  2 
+  45  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 
+  60  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 
+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  sin(wk) 

2  2  2 
+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  sin(wk) 

2  2 
+  180  cos  (apk)  sin(apk)  cos(fk)  cos  (jk)  sin(wk) 

3  2 

-  18  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(vk) 

2 

+  18  cos(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(wk) 

3 

-  90  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

+  90  cos(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

3  2 

-  48  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

2 

+  48  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  'sin(wk) 

3 

-  120  cos  (apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

+  120  cos(apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

3 

-  90  cos  (apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(wk) 

+  90  cos(apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(wk) 

3  2 

+  66  cos  (apk)  ek  sin(fk)  cos(jk)  sin(wk) 

2 

-  66  cos(apk)  ek  sin(fk)  cos(jk)  sin(wk) 

3 

+  120  cos  (apk)  sin(fk)  cos(jk)  sin(vk)  -  120  cos(apk)  sin(fk)  cos(jk)  sin(sk) 
2  2 

+  3  cos  (apk)  sin(apk)  ek  cos(5  fk)  sin(wk) 

2  2 

-  3  sin(apk)  ek  cos(5  fk)  sin(sk)  +  15  cos  (apk)  sin(apk)  ek  cos(4  fk) 

sin(vk)  -  15  sin(apk)  ek  cos(4  fk)  sin(sk) 

2  2 

-  5  cos  (apk)  sin(apk)  ek  cos(3  fk)  sin(wk) 

2  2 

+  5  sin(apk)  ek  cos(3  fk)  sin(wk)  +  20  cos  (apk)  sin(apk)  cos(3  fk)  sin(wk) 

2 

-  20  sin(apk)  cos(3  fk)  sin(wk)  -  60  cos  (apk)  sin(apk)  ek  cos(2  fk)  sin(wk) 
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2  2 

+  60  sin(apk)  ek  cos(2  fk)  sin(wk)  -  30  cos  (apk)  sin(apk)  ek  cos(fk)  sin(wk) 
2  2 

+  30  sin(apk)  ek  cos(fk)  sin(wk)  -  180  cos  (apk)  sin(apk)  cos(fk)  sin(wk) 

3  2  3  3 

+  180  sin(apk)  cos(fk)  sin(wk)  -  12  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk) 

3  3  3 

-  60  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk) 

3  2  3  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3  3 

-  80  cos  (apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3  3 

-  60  cos  (apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  2  2  3 

-  72  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  360  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  2  2  3 

-  312  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  480  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  1080  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

2  2.23 

-  336  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  240  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos  (wk) 

3  2  3 

-  36  cos  (apk)  ek  cos(5  fk)  cos(jk)  cos  (wk) 

2  3 

+  36  cos(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk) 

3  3 

-  180  cos  (apk)  ek  cos(4  fk)  cos(jk)  cos  (wk) 

3 

+  180  cos(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk) 

3  2  3 

-  120  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

+  120  cos(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 

3  3 

-  240  cos  (apk)  cos(3  fk)  cos(jk)  cos  (wk) 

3 

+  240  cos(apk)  cos(3  fk)  cos(jk)  cos  (wk) 

3  3 

-  360  cos  (apk)  ek  cos(2  fk)  cos(jk)  cos  (wk) 

3 

+  360  cos(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk) 

3  2  3 

-  180  cos  (apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

2  3 

+  180  cos(apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

2  2  3 

-  24  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (wk) 

2  3 

+  24  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (wk) 
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2  3 

-  120  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (wk) 

3 

+  120  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (wk) 

2  2  3 

-  104  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (wk) 

2  3 

+  104  sin(apk)  ek  sm(fk)  cos(2  fk)  cos  (wk) 

2  3 

-  160  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (wk) 

3 

+  160  sin(apk)  sin(fk)  cos(2  fk)  cos  (wk) 

2  3 

-  360  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (wk) 

3 

+  360  sin(apk)  ek  cos(fk)  sin(fk)  cos  (wk) 

2  2  3 

-  112  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (wk) 

2  3  2  3 

+  112  sin(apk)  ek  sin(fk)  cos  (wk)  -  80  cos  (apk)  sin(apk)  sin(fk)  cos  (wk) 

3  3  2  3 

+  80  sin(apk)  sin(fk)  cos  (wk)  +  9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos(wk) 

3  3 

+  45  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

3  2  3 

+  45  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

3  3 

+  60  cos  (apk)  cos(3  fk)  cos  (jk)  cos(wk) 

3  3 

+  180  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

3  2  3 

+  90  cos  (apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

3  3 

+  1 80  cos  (apk)  cos(fk)  cos  (jk)  cos(wk) 

2  2  2 
+  54  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

2  2 
+  270  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  2  2 

+  264  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 
+  360  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 
+  990  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(wk) 

2  2  2 
+  402  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos(wk) 

2  2 
+  360  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos(wk) 

3  2 

+  27  cos  (apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

2 

-  27  cos(apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

3 

+  135  cos  (apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

-  135  cos(apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

3  2 

+  75  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 
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2 

-  75  cos(apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

3 

+  180  cos  (apk)  cos(3  fk)  cos(jk)  cos(wk) 

-  180  cos(apk)  cos(3  fk)  cos(jk)  cos(wk) 

3 

+  180  cos  (apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

-  180  cos(apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

3  2 

+  90  cos  (apk)  ek  cos(fk)  cos(jk)  cos(wk) 

2 

-  90  cos(apk)  ek  cos(fk)  cos(jk)  cos(wk) 

3 

-  180  cos  (apk)  cos(fk)  cos(jk)  cos(wk)  +  180  cos(apk)  cos(fk)  cos(jk)  cos(wk) 

2  2 

+  18  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(wk) 

2 

-  18  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(wk) 

2 

+  90  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(wk) 

-  90  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(wk) 

2  2 

+  48  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(wk) 

2 

-  48  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(wk) 

2 

+  120  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(wk) 

-  120  sin(apk)  sin(fk)  cos(2  fk)  cos(wk) 

2 

+  90  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(wk) 

-  90  sin(apk)  ek  cos(fk)  sin(fk)  cos(wk) 

2  2  2 

-  66  cos  (apk)  sin(apk)  ek  sin(fk)  cos(wk)  +  66  sin(apk)  ek  sin(fk)  cos(wk) 

2 

-  120  cos  (apk)  sin(apk)  sin(fk)  cos(wk)  +  120  sin(apk)  sin(fk)  cos(wk))  zb2 
/ (240  nk) 


^(ra25l522)  = 

2  2  3  2 

n  (24  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 
2  3  2 

+  120  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  3  2 

+  104  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  160  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  360  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 
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2  2  3  2 

+  112  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  3  2 

+  80  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  36  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  24  cos(apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  180  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  120  cos(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  120  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  80  cos(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk)  sin(vk) 

3  2  2 

♦  240  cos  (apk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  160  cos(apk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2 

+  360  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  240  cos(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

3  2  2  2 

+  180  cos  (apk)  ek  cos(fk)  cos  (jk)  cos  (wk)  sin(vk) 

2  2  2 

-  120  cos(apk)  ek  cos(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  72  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  24  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  360  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk)  sizi(wk) 

2 

-  120  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  312  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 
2  2 

-  104  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  480  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

*  160  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  ’  2 
+  1080  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  360  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  336  cos  (apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  112  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  240  cos  (apk)  sin(apk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  80  sin(apk)  sin(fk)  cos(jk)  cos  (wk)  sin(wk) 

3  2  2 

+  12  cos  (apk)  ek  cos(5  fk)  cos  (wk)  sin(wk) 
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2  2 

-  12  cos(apk)  ek  cos(5  fk)  cos  (sk)  sin(ffk) 

3  2 

+  60  cos  (apk)  ek  cos(4  fk)  cos  (wk)  sin(wk) 

2 

-  60  cos(apk)  ek  cos(4  fk)  cos  (wk)  sin(wk) 

3  2  2 

+  40  cos  (apk)  ek  cos(3  fk)  cos  (sk)  sin(wk) 

2  2 

-  40  cos(apk)  ek  cos(3  fk)  cos  (sk)  sin(wk) 

3  2 

+  80  cos  (apk)  cos(3  fk)  cos  (wk)  sin(wk) 

2 

-  80  cos(apk)  cos(3  fk)  cos  (wk)  sin(sk) 

3  2 

+  120  cos  (apk)  ek  cos(2  fk)  cos  (sk)  sin(wk) 

2 

-  120  cos(apk)  ek  cos(2  fk)  cos  (wk)  sin(wk) 

3  2  2 

+  60  cos  (apk)  ek  cos(fk)  cos  (wk)  sin(wk) 

2  2 

-  60  cos(apk)  ek  cos(fk)  cos  (sk)  sin(wk) 

2  2  3 

-  6  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 

2  3 

-  30  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2  3 

-  56  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  3 

-  40  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(Hk) 

2  3 

-  270  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  3 

-  178  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2  3 

-  200  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  sin(wk) 

3  2  2 

-  9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

2  2 

+  6  cos(apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

3  2 

-  45  cos  (apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

2 

+  30  cos(apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

3  2  2 

-  45  cos  (apk)  ek  •  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 

+  30  cos(apk)  ek  cos(3  fk)  cos  (jk)  sin(wk) 

3  2 

-  60  cos  (apk)  cos (3  fk)  cos  (jk)  sin(wk) 

2 

+  40  cos(apk)  cos(3  fk)  cos  (jk)  sin(wk) 

3  2 

-  180  cos  (apk)  ek  cos(2  fk)  cos  (jk)  sin(wk) 

2 

+  120  cos(apk)  ek  cos(2  fk)  cos  (jk)  sin(wk) 

3  2  2 

-  90  cos  (apk)  ek  cos(fk)  cos  (jk)  sin(wk) 
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2  2 

+  60  cos(apk)  ek  cos(fk)  cos  (jk)  sin(wk) 

3  2 

-  180  cos  (apk)  cos(fk)  cos  (jk)  sin(wk) 

2 

+  120  cos(apk)  cos(fk)  cos  (jk)  sin(wk) 

2  2 

-  18  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(wk) 

2 

+  6  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(sk) 

2 

-  90  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

+  30  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

2  2 

-  48  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

2 

+  16  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(vk) 

2 

-  120  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

*  40  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(sk) 

2 

-  90  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(yk) 

+  30  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(yk) 

2  2 

+  66  cos  (apk)  sin(apk)  ek  sin(fk)  cos(jk)  sm(wk) 

2 

-  22  sin(apk)  ek  sin(fk)  cos(jk)  sin(vk) 

2 

+  120  cos  (apk)  sin(apk)  sin(fk)  cos(jk)  sin(wk) 

3  2 

-  40  sin(apk)  sm(fk)  cos(jk)  sin(sk)  -  3  cos  (apk)  ek  cos(5  fk)  sin(vk) 

2  3 

+  3  cos(apk)  ek  cos(5  fk)  sin(yk)  -  15  cos  (apk)  ek  cos(4  fk)  sin(sk) 

3  2 

+  15  cos(apk)  ek  cos(4  fk)  sin(wk)  +  5  cos  (apk)  ek  cos(3  fk)  sin(wk) 

2  3 

-  5  cos(apk)  ek  cos(3  fk)  sin(wk)  -  20  cos  (apk)  cos(3  fk)  sin(sk) 

3 

+  20  cos(apk)  cos(3  fk)  sin(vk)  +  60  cos  (apk)  ek  cos(2  fk)  sin(sk) 

3  2 

-  60  cos(apk)  ek  cos(2  fk)  sin(wk)  +  30  cos  (apk)  ek  cos(fk)  sin(sk) 

2  3 

-  30  cos(apk)  ek  cos(fk)  sin(wk)  +  180  cos  (apk)  cos(fk)  sin(wk) 

2  2  3 

-  180  cos(apk)  cos(fk)  sin(wk)  -  12  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk) 

3  2  3  3 

cos  (wk)  -  60  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (*k) 

2  2  3  3 

-  40  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

-  80  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

-  120  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (sk) 

2  2  3  3 

-  60  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (yk) 
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3  2  2  3 

72  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3 

48  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

3  2  3 

360  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3 

240  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  2  2  3 

312  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 
2  2  3 

208  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3 

480  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  3 

320  cos(apk)  sm(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3 

1080  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

2  3 

720  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 

3  2  2  3 

336  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (wk) 

2  2  3 

224  cos(apk)  ek  sin(fk)  cos  (jk)  cos  (wk) 

3  2  3 

240  cos  (apk)  sin(fk)  cos  (jk)  cos  (wk) 

2  3 

160  cos(apk)  sin(fk)  cos  (jk)  cos  (wk) 

2  2  3 

36  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk) 

2  3 

12  sin(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk) 

2  3 

180  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk) 

3 

60  sin(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk) 

2  2  3 
120  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 
2  3 

40  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

240  cos  (apk)  sin(apk)  cos(3  fk)  cos(jk)  cos  (wk) 

3 

80  sin(apk)  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 
360  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk) 

3 

120  sin(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk) 

2  2  3 
180  cos  (apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

2  3 

60  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk) 

3  2  3 

24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (wk) 

2  3 

24  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (wk) 

3  3 
120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (wk) 


58 


3 

-  120  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (wk) 

3  2  3 

+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (wk) 

2  3 

-  104  cos(apk)  ek  sirt(fk)  cos(2  fk)  cos  (wk) 

3  3 

+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (wk) 

3 

-  160  cos(apk)  sin(fk)  cos(2  fk)  cos  (wk) 

3  3 

+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (wk) 

3 

-  360  cos(apk)  ek  cos(fk)  sin(fk)  cos  (wk) 

3  2  3  2  3 

+  112  cos  (apk)  ek  sin(fk)  cos  (wk)  -  112  cos(apk)  ek  sin(fk)  cos  (wk) 
3  3  3 

+  80  cos  (apk)  sin(fk)  cos  (wk)  -  80  cos(apk)  sin(fk)  cos  (wk) 

2  2  3 

+  9  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos(wk) 

2  3 

+  45  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

2  2  3 

+  45  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

2  3 

+  60  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  3 

+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

2  2  3 

+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

2  3 

+  180  cos  (apk)  sin(apk)  cos(fk)  cos  (jk)  cos(wk) 

3  2  2 

-  54  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

2  2 

+  36  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

3  2 

-  270  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(wk) 

2 

+  180  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  co6(wk) 

3  2  2 

-  264  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 

+  176  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

3  2 

-  360  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2 

+  240  cos(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

3  2 

-  990  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(wk) 

2 

+  660  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(wk) 

3  2  2 

-  402  cos  (apk)  ek  sin(fk)  cos  (jk)  cos(wk) 

2  2 

+  268  cos(apk)  ek  sin(fk)  cos  (jk)  cos(wk) 

3  2 

-  360  cos  (apk)  sin(fk)  cos  (jk)  cos(wk) 
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2 

+  240  cos(apk)  sin(fk)  cos  (jk)  cos(wk) 

2  2 

+  27  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

2 

-  9  sin(apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

2 

+  135  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

-  45  sin(apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

2  2 

+  75  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos(sk) 

2 

-  25  sin(apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  180  cos  (apk)  sin(apk)  cos(3  fk)  cos(jk)  cos(wk) 

-  60  sin(apk)  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  cos(sk) 

-  60  sin(apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

2  2 

+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos(wk) 

2 

-  30  sin(apk)  ek  cos(fk)  cos(jk)  cos(wk) 

2 

-  180  cos  (apk)  sin(apk)  cos(fk)  cos(jk)  cos(wk) 

3  2 

+  60  sin(apk)  cos(fk)  cos(jk)  cos(wk)  -  18  cos  (apk)  ek  sin(fk)  cos(4  fk) 

2 

cos(wk)  +  18  cos(apk)  ek  sin(fk)  cos(4  fk)  cos(wk) 

3 

-  90  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(wk) 

+  90  cos(apk)  ek  sin(fk)  cos(3  fk)  cos(vk) 

3  2 

-  48  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(wk) 

2 

+  48  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(sk) 

3 

-  120  cos  (apk)  sin(fk)  cos(2  fk)  cos(wk) 

+  120  cos(apk)  sin(fk)  cos(2  fk)  cos(wk) 

3 

-  90  cos  (apk)  ek  cos(fk)  sin(fk)  cos(sk) 

3  2 

+  90  cos(apk)  ek  cos(fk)  sin(fk)  cos(wk)  +  66  cos  (apk)  ek  sin(fk)  cos(wk) 
2  3 

-  66  cos(apk)  ek  sin(fk)  cos(sk)  +  120  cos  (apk)  sin(fk)  cos(wk) 

-  120  cos(apk)  sin(fk)  cos(sk))  zb2/(240  nk) 
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6(za2S2Sl2)  = 

3  2  3  2 

-  n  (24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (vk)  sin(vk) 
2  3  2 

-  24  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  3  2 

+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2 

-  120  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2  3  2 

+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  3  2 

-  104  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  3  2 

+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2 

-  160  cos(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  3  2 

+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2 

-  360  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2  3  2 

+  112  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (vk)  sin(wk) 

2  3  2 

-  112  cos(apk)  ek  sin(fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  3  2 

+  80  cos  (apk)  sin(fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2 

-  80  cos(apk)  sin(fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2  2 

-  36  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

+  12  sin(apk)  ek  cos(S  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

-  180  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2 

+  60  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2  2 

-  120  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

+  40  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

-  240  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2 

+  80  sin(apk)  cos(3  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

-  360  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2 

+  120  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2  2 

-  180  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (vk)  sin(vk) 

2  2  2 

+  60  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (vk)  sin(vk) 

3  2  2 

+  72  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (vk)  sin(vk) 

2  2 

-  48  cos(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (vk)  sin(vk) 
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3  2 

+  360  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (yk)  sin(yk) 

2 

-  240  cos(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (yk)  sin(yk) 

3  2  2 

+  312  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (yk)  sin(yk) 
2  2 

-  208  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (yk)  sin(yk) 

3  2 

+  480  cos  (apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (yk)  sin(yk) 

2 

-  320  cos(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (yk)  sin(yk) 

3  2 

+  1080  cos  (apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

2 

-  720  cos(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

3  2  2 

+  336  cos  (apk)  ek  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

2  2 

-  224  cos(apk)  ek  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

3  2 

+  240  cos  (apk)  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

2 

-  160  cos(apk)  sin(fk)  cos(jk)  cos  (yk)  sin(yk) 

2  2  2 

-  12  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (yk)  sin(yk) 

2  2 

-  60  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (yk)  sin(yk) 

2  2  2 

-  40  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (yk)  sin(yk) 

2  2 

-  80  cos  (apk)  sin(apk)  cos(3  fk)  cos  (yk)  sin(yk) 

2  2 

-  120  cos  (apk)  sin(apk)  ek  cos  (2  fk)  cos  (yk)  sin(wk) 

2  2  2 

-  60  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (yk)  sin(yk) 

3  2  3 

-  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(yk) 

2  3 

+  6  cos(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(yk) 

3  3 

-  30  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(yk) 

3 

+  30  cos(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(yk) 

3  2  3 

-  56  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(yk) 

2  3 

+  56  cos(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(yk) 

3  3 

-  40  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(yk) 

3 

+  40  cos(apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(yk) 

3  3 

-  270  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(yk) 

3 

+  270  cos(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(yk) 

3  2  3 

-  178  cos  (apk)  ek  sin(fk)  cos  (jk)  sin(yk) 
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2  3 

+  178  cos(apk)  ek  sin(fk)  cos  (jk)  sin(sk) 

3  3 

-  200  cos  (apk)  sin(fk)  cos  (jk)  sin(sk) 

3 

+  200  cos(apk)  sirv(fk)  cos  (jk)  sin(wk) 

2  2  2 
+  9  cos  (apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  sin(trk) 

2  2 

-  3  sin(apk)  ek  cos(5  fk)  cos  (jk)  sin(wk) 

2  2 
+  45  cos  (apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

2 

-  15  sin(apk)  ek  cos(4  fk)  cos  (jk)  sin(wk) 

2  2  2 
+  45  cos  (apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  sin(sk) 
2  2 

-  15  sin(apk)  ek  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 
*►  60  cos  (apk)  sin(apk)  cos(3  fk)  cos  (jk)  sin(sk) 

2 

-  20  sin(apk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2 
+  180  cos  (apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  sin(wk) 

2 

-  60  sin(apk)  ek  cos(2  fk)  cos  (jk)  sin(vk) 

2  2  2 
+  90  cos  (apk)  sin(apk)  ek  cos(fk)  cos  (jk)  sin(sk) 

2  2 

-  30  sm(apk)  ek  cos(fk)  cos  (jk)  sin(wk) 

2  2 
*  180  cos  (apk)  sin(apk)  cos(fk)  cos  (jk)  sin(wk) 

2 

-  60  sm(apk)  cos(fk)  cos  (jk)  sin(*k) 

3  2 

-  18  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(wk) 

2 

+  12  cos(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  sin(wk) 

3 

-  90  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

+  60  cos(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  sin(wk) 

3  2 

-  48  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

2 

+  32  cos(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

3 

-  120  cos  (apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

+  80  cos(apk)  sin(fk)  cos(2  fk)  cos(jk)  sin(wk) 

3 

-  90  cos  (apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(wk) 

+  60  cos(apk)  ek  cos(fk)  sin(fk)  cos(jk)  sin(wk) 

3  2 

+  66  cos  (apk)  ek  sin(fk)  cos(jk)  sin(wk) 

2 

-  44  cos(apk)  ek  sin(fk)  cos(jk)  sin(wk) 
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3 

+  120  cos  (apk)  sin(fk)  cos(jk)  sin(wk)  -  80  cos(apk)  sin(fk)  cos(jk)  sin(wk) 
2  2 

+  3  cos  (apk)  sin(apk)  ek  cos(5  fk)  sin(sk) 

2 

+  15  cos  (apk)  6in(apk)  ek  cos(4  Ik)  sin(wk) 

2  2 

-  5  cos  (apk)  sin(apk)  ek  cos(3  fk)  sin(wk) 

2 

+  20  cos  (apk)  sin(apk)  cos(3  fk)  sin(wk) 

2 

-  60  cos  (apk)  sin(apk)  ek  cos(2  fk)  sin(wk) 

2  2 

-  30  cos  (apk)  sin(apk)  ek  cos(fk)  sin(wk) 

2 

-  180  cos  (apk)  sin(apk)  cos(fk)  6in(wk) 

3  2  3  3 

-  12  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk) 

2  3  3 

+  12  cos(apk)  ek  cos(5  fk)  co6  (jk)  cos  (wk) 

3  3  3 

-  60  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk) 

3  3 

+  60  cos(apk)  ek  cos(4  fk)  cos  (jk)  cos  (wk) 

3  2  3  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3  3 

+  40  cos(apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3  3 

-  80  cos  (apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3 

+  80  cos(apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

3  3  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

3  3 

+  120  cos(apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

3  2  3  3 

-  60  cos  (apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  3  3 

+  60  cos(apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  2  2  3 

-  72  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3 

+  24  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  360  cos  (apk)  sin(apk)  ek  6in(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  3 

+  120  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  2  2  3 

-  312  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

+  104  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  480  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  3 

+  160  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  1080  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk) 
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2  3 

360  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (sk) 

2  2  2  3 

336  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (sk) 

2  2  3 

112  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (sk) 

2  2  3 

240  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos  (sk) 

2  3 

80  sin(apk)  sin(fk)  cos  (jk)  cos  (sk) 

3  2  3 

36  cos  (apk)  ek  cos(5  fk)  cos(jk)  cos  (sk) 

2  3 

24  cos(apk)  ek  cos(5  fk)  cos(jk)  cos  (sk) 

3  3 
180  cos  (apk)  ek  cos(4  fk)  cos(jk)  cos  (sk) 

3 

120  cos(apk)  ek  cos(4  fk)  cos(jk)  cos  (sk) 

3  2  3 

120  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos  (sk) 

2  3 

80  cos(apk)  ek  cos(3  fk)  cos(jk)  cos  (sk) 

3  3 

240  cos  (apk)  cos(3  fk)  cos(jk)  cos  (sk) 

3 

160  cos(apk)  cos(3  fk)  cos(jk)  cos  (sk) 

3  3 

360  cos  (apk)  ek  cos(2  fk)  cos(jk)  cos  (sk) 

3 

240  cos(apk)  ek  cos(2  fk)  cos(jk)  cos  (sk) 

3  2  3 

180  cos  (apk)  ek  cos(fk)  cos(jk)  cos  (sk) 

2  3 

120  cos(apk)  ek  cos(fk)  cos(jk)  cos  (sk) 

2  2  3 

24  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (sk) 

‘  2  3 

120  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (sk) 

2  2  3 

104  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (sk) 
2  3 

160  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (sk) 

2  3 

360  cos  (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (sk) 

2  2  3 

112  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (sk) 

2  3 

80  cos  (apk)  sin(apk)  sin(fk)  cos  (sk) 

3  2  3 

9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos(sk) 

2  3 

9  cos(apk)  ek  cos(5  fk)  cos  (jk)  cos(sk) 

3  3 

45  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos(sk) 

3 

45  cos(apk)  ek  cos(4  fk)  cos  (jk)  cos(sk) 

3  2  3 

45  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(sk) 
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2  3 

-  45  cos(apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

3  3 

+  60  cos  (apk)  cos(3  fk)  cos  (jk)  cos(wk) 

3 

-  60  cos(apk)  cos(3  fk)  cos  (jk)  cos(wk) 

3  3 

+  180  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos(sk) 

3 

-  180  cos(apk)  ek  cos(2  fk)  cos  (jk)  cos(sk) 

3  2  3 

+  90  cos  (apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

2  3 

-  90  cos(apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

3  3 

+  180  cos  (apk)  cos(fk)  cos  (jk)  cos(wk) 

3 

-  180  cos(apk)  cos(fk)  cos  (jk)  cos(wk) 

2  2  2 
+  54  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

2  2 

-  18  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos(wk) 

2  2 
+  270  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(vk) 

2 

-  90  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  2  2 
+  264  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 
2  2 

-  88  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 
+  360  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2 

-  120  6in(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 
♦  990  cos  (apk)  sm(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(yk) 

2 

-  330  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos(sk) 

2  2  2 
+  402  cos  (apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos(yk) 

2  2 

-  134  sin(apk)  ek  sin(fk)  cos  (jk)  cos(yk) 

2  2 
+  360  cos  (apk)  sin(apk)  sin(fk)  cos  (jk)  cos(wk) 

2 

-  120  sin(apk)  sin(fk)  cos  (jk)  cos(wk) 

3  2 

+  27  cos  (apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

2 

-  18  cos(apk)  ek  cos(5  fk)  cos(jk)  cos(wk) 

3 

+  135  cos  (apk)  ek  cos(4  fk)  cos(jk)  co6(wk) 

-  90  cos(apk)  ek  cos(4  fk)  cos(jk)  cos(wk) 

3  2 

+  75  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos(wk) 

2 

-  50  cos(apk)  ek  cos(3  fk)  cos(jk)  cos(uk) 
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3 

+  180  cos  (apk)  cos(3  fk)  cos(jk)  cos(wk) 

-  120  cos(apk)  cos(3  fk)  cos(jk)  cos(wk) 

3 

+  180  cos  (apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

-  120  cos(apk)  ek  cos(2  fk)  cos(jk)  cos(wk) 

3  2 

+  90  cos  (apk)  ek  cos(fk)  cos(jk)  cos(wk) 

2 

-  60  cos(apk)  ek  cos(fk)  cos(jk)  cos(wk) 

3 

-  180  cos  (apk)  cos(fk)  cos(jk)  cos(wk)  +  120  cos(apk)  cos(fk)  cos(jk)  cos(wk) 

2  2 

+  18  cos  (apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(wk) 

2 

+  90  cos  (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(wk) 

2  2 

+  48  cos  (apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(wk) 

2 

+  120  cos  (apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(wk) 

2 

+  90  cos  (apk)  sm(apk)  ek  cos(fk)  sin(fk)  cos(vk) 

2  2 

-  66  cos  (apk)  sin(apk)  ek  sin(fk)  cos(wk) 

2 

-  120  cos  (apk)  sin(apk)  sin(fk)  cos(wk))  zb2/(240  nk) 


£(ra2SlS2S3)  = 

2  2  2 

-  sin(jk)  n  (24  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk) 

2  2 

sin(wk)  +  120  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

+  104  cos(apk)  sm(apk)  ek  sm(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

+  160  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

+  360  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  siii(wk) 

2  2  2 

+  112  cos(apk)  sin(apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

+  80  cos(apk)  sin(apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

+  24  cos  (apk)  ek  cos(5  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  12  ek  cos(5  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

+  120  cos  (apk)  ek  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  60  ek  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 

+  80  cos  (apk)  ek  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 
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2  2 

-  40  ek  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  160  cos  (apk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  80  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 
+  240  cos  (apk)  ek  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  120  ek  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  120  cos  (apk)  ek  cos(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  60  ek  cos(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

+  24  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (wk)  sin(wk) 

2 

+  120  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos  (wk)  sin(wk) 

2  2 

+  104  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2 

+  160  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (wk)  sxn(wk) 

2 

+  360  cos (apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (wk)  sin(wk) 

2  2 

+  112  cos(apk)  sm(apk)  ek  sin(fk)  cos  (wk)  sin(wk) 

2 

+  80  cos (apk)  sin(apk)  sin(fk)  cos  (wk)  sin(wk) 

2  2 

-  6  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 

2 

-  30  cos(apk)  sin(apk)  ek  sin(fk)  cos (3  fk)  cos  (jk)  sin(wk) 

2  2 

-  56  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2 

-  40  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2 

-  270  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2  2 

-  178  cos(apk)  sin(apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2 

-  200  cos(apk)  sin(apk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  2 

-  6  cos  (apk)  ek  cos(5  fk)  cos(jk)  sin(wk)  +  3  ek  cos(5  fk)  cos(jk)  sin(wk) 

2 

-  30  cos  (apk)  ek  cos(4  fk)  cos(jk)  sin(wk)  +  15  ek  cos(4  fk)  cos(jk)  sin(wk) 

2  2 

-  30  cos  (apk)  ek  cos(3  fk)  cos(jk)  sin(wk) 

2  2 

+  15  ek  cos(3  fk)  cos(jk)  sin(wk)  -  40  cos  (apk)  cos(3  fk)  cos(jk)  sin(wk) 

2 

+  20  cos (3  fk)  cos(jk)  sin(wk)  -  120  cos  (apk)  ek  cos(2  fk)  cos(jk)  sin(wk) 

2  2 

+  60  ek  cos(2  fk)  cos(jk)  sin(wk)  -  60  cos  (apk)  ek  cos(fk)  cos(jk)  sin(wk) 

2  2 

+  30  ek  cos(fk)  cos(jk)  sin(wk)  -  120  cos  (apk)  cos(fk)  cos(jk)  sin(wk) 

2 

+  60  cos(fk)  cos(jk)  sin(wk)  -  6  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk) 
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d 


sin(sk)  -  30  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  sin(sk) 

2 

-  16  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  sin(sk) 

-  40  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  sin(sk) 

-  30  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  sin(sk) 

2 

+  22  cos(apk)  sin(apk)  ek  Bin(fk)  sin(sk) 


+  40  cos(apk)  sin(apk)  sin(fk)  sin(sk) 

2  2  3 

-  12  cos(apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos  (sk) 

2  3 

-  60  cos(apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cob  (sk) 

2  2  3 

-  40  cos(apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos  (sk) 

2  3 

-  80  cos(apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos  (sk) 

2  3 

-  120  cos(apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos  (sk) 

2  2  3 

-  60  cos(apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos  (sk) 

2  2  3 

+  48  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (sk) 

2  3 

-  24  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (sk) 

2  3 

+  240  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (sk) 

3 

-  120  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (sk) 

2  2  3 

+  208  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (sk) 
2  3 

-  104  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (sk) 

2  3 

+  320  cos  (apk)  sin(fk)  cob(2  fk)  cos(jk)  cos  (sk) 

3 

-  160  sin(fk)  cos(2  fk)  cos(jk)  cos  (sk) 

2  3 


+  720  cos  (apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (sk) 

3 

-  360  ek  cos(fk)  sin(fk)  cos(jk)  cos  (sk) 

2  2  3 

+  224  cos  (apk)  ek  sin(fk)  cos(jk)  cos  (sk) 

2  3  2  3 

-  112  ek  sin(fk)  cos(jk)  cos  (sk)  +  160  cos  (apk)  sin(fk)  cos(jk)  cos  (sk) 

3  2  3 

-  80  sin(fk)  cos(jk)  cos  (sk)  -  12  cos(apk)  sin(apk)  ek  cos(5  fk)  cos  (sk) 

3 

-  60  cos(apk)  sin(apk)  ek  cos(4  fk)  cos  (sk) 

2  3 

-  40  cos(apk)  sin(apk)  ek  cos(3  fk)  cos  (sk) 

3 

-  80  cos(apk)  sin(apk)  cos(3  fk)  cos  (sk) 

3 

-  120  cos(apk)  sin(apk)  ek  cos(2  fk)  cos  (sk) 

2  3 

-  60  cos(apk)  sin(apk)  ek  cos(fk)  cos  (sk) 
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2  2 

+  9  cos(apk)  sin(apk)  ek  cos(5  fk)  cos  (jk)  cos(yk) 

2 

+  45  cos(apk)  sin(apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

2  2 

+  45  cos(apk)  sin(apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

2 

+  60  cos(apk)  sin(apk)  cos(3  fk)  cos  (jk)  cos(wk) 

2 

+  180  cos(apk)  sin(apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

2  2 

+  90  cos(apk)  sin(apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

2 

+  180  cos(apk)  sin(apk)  cos(fk)  cos  (jk)  cos(wk) 

2  2 

-  36  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos(sk) 

2 

+  18  ek  sin(fk)  cos(4  fk)  cos(jk)  cos(wk) 

2 

-  180  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos(wk) 

+  90  ek  sin(fk)  cos(3  fk)  cos(jk)  cos(wk) 

2  2 

-  176  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

2 

+  88  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

2 

-  240  cos  (apk)  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

+  120  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

2 

-  660  cos  ’(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos(wk) 

+  330  ek  cos(fk)  sin(fk)  cos(jk)  cos(wk) 

2  2  2 

-  268  cos  (apk)  ek  sin(fk)  cos(jk)  cos(wk)  +  134  ek  sin(fk)  cos(jk)  cos(yk) 

2 

-  240  cos  (apk)  sin(fk)  cos(jk)  cos(wk)  +  120  sin(fk)  cos(jk)  cos(wk) 

2 

+  9  cos(apk)  sin(apk)  ek  cos(5  fk)  cos(wk) 


+  45  cos(apk)  sin(apk)  ek  cos(4  fk)  cos(wk) 

2 

+  25  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(wk) 

+  60  cos(apk)  sin(apk)  cos(3  fk)  cos(wk) 

+  60  cos(apk)  sin(apk)  ek  cos(2  fk)  cos(wk) 

2 

+  30  cos(apk)  sin(apk)  ek  cos(fk)  cos(wk) 

-  60  cos(apk)  sin(apk)  cos(fk)  cos(wk))  zb2/(240  nk) 
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6(:a2S3S22)  = 


2  2  2  2 
-  sin(jk)  n  (24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 
2  2  2 
+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 
+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  360  cos  (apk)  ck  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 
+  112  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  80  cos  (apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  24  cos(apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  120  cos(apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  80  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  160  cos(apk)  sin(apk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  240  cos(apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk)  sm(wk) 

2  2 

-  120  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (wk)  sin(wk) 

2  2 

-  24  ek  sin(fk)  cos(4  fk)  cos  (wk)  sin(wk) 

2  2 
+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (wk)  sm(wk) 

2 

-  120  ek  sin(fk)  cos(3  fk)  cos  (wk)  sin(wk) 

2  2  2 
+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2  2 

-  104  ek  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2  2 
+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2 

-  160  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2  2 
+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (wk)  sin(wk) 

2 

-  360  ek  cos(fk)  sin(fk)  cos  (wk)  sin(wk) 

2  2  2 
+  112  cos  (apk)  ek  sin(fk)  cos  (wk)  sin(wk) 

2  2  2  2 
*  112  ek  sin(fk)  cos  (wk)  sin(wk)  *  80  cos  (apk)  sin(fk)  cos  (wk)  sin(wk) 

2  2  2  2 

-  80  sin(fk)  cos  (wk)  sin(wk)  -  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk) 

2  2 
sin(wk)  -  30  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2  2 

-  56  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 
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2  2 

-  40  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  2 

-  270  cos  (apk)  ck  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  2 

-  178  cos  (apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2  2 

-  200  cos  (apk)  sin(fk)  cos  (jk)  sin(wk) 

2 

+  6  cos(apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  sin(wk) 

+  30  cos (apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  sin(wk) 

2 

+  30  cos(apk)  sin(apk)  ek  cos (3  fk)  cos(jk)  sin(wk) 

+  40  cos(apk)  sin(apk)  cos(3  fk)  cos(jk)  sin(wk) 

+  120  cos(apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  sin(wk) 

2 

+  60  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  sin(wk) 

+  120  cos(apk)  sin(apk)  cos(fk)  cos(jk)  sin(wk) 

2  2  2 

-  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  sin(wk)  +  6  ek  sin(fk)  cos(4  fk)  sin(wk) 

2 

-  30  cos  (apk)  ek  sin(fk)  cos(3  fk)  sin(wk)  +  30  ek  sin(fk)  cos(3  fk)  sin(wk) 

2  2 

-  16  cos  (apk)  ek  sin(fk)  cos(2  fk)  sin(wk) 

2  2 

+  16  ek  sin(fk)  cos(2  fk)  sin(wk)  -  40  cos  (apk)  sin(fk)  cos(2  fk)  sin(wk) 

2 

+  40  sin(fk)  cos(2  fk)  sin(wk)  -  30  cos  (apk)  ek  cos(fk)  sin(fk)  sin(wk) 

2  2 

+  30  ek  cos(fk)  sin(fk)  sin(wk)  +  22  cos  (apk)  ek  sin(fk)  sin(wk) 

2  2 

-  22  ek  sin(fk)  sin(wk)  +  40  cos  (apk)  sin(fk)  sin(wk)  -  40  sin(fk)  sin(yk) 

2  2  2  3 

-  12  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  60  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos  (sk) 

2  2  2  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  80  cos  (apk)  cos(3  fk)  cos  (jk)  cos  (wk) 

2  2  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

2  2  2  3 

-  60  cos  (apk)  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  3 

-  48  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk) 

3 

-  240  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

-  208  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 

3 

-  320  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 

3 

-  720  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk) 
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2  3 

-  224  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk) 

3 

-  160  cos(apk)  sin(apk)  sin(fk)  cos(jk)  cos  (wk) 

2  2  3  2  3 

-  12  cos  (apk)  ek  cos(5  fk)  cos  (wk)  +  12  ek  cos(S  fk)  cos  (wk) 

2  3  3 

-  60  cos  (apk)  ek  cos(4  fk)  cos  (wk)  +  60  ek  cos(4  fk)  cos  (wk) 

2  2  3  2  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (wk)  +  40  ek  cos(3  fk)  cos  (wk) 

2  3  3 

-  80  cos  (apk)  cos(3  fk)  cos  (wk)  +  80  cos(3  fk)  cos  (wk) 

2  3  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (wk)  +  120  ek  cos(2  fk)  cos  (wk) 

2  2  3  2  3 

-  60  cos  (apk)  ek  cos(fk)  cos  (wk)  +  60  ek  cos(fk)  cos  (wk) 

2  2  2 
♦  9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos(wk) 

2  2 
+  45  cos  (apk)  ek  cos(4  fk)  cos  (jk)  cos(wk) 

2  2  2 
+  45  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

2  2 
+  60  cos  (apk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  2 
+  180  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

2  2  2 
+  90  cos  (apk)  ek  cos(fk)  cos  (jk)  cos(wk) 

2  2 

+  180  cos  (apk)  cos(fk)  cos  (jk)  cos(wk) 

2 

+  36  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos(wk) 

+  180  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  176  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

+  240  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

+  660  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos(wk) 

2 

+  268  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos(wk) 

+  240  cos(apk)  sin(apk)  sm(fk)  cos(jk)  cos(wk) 

2  2  2 
+  9  cos  (apk)  ek  cos(5  fk)  cos(wk)  -  9  ek  cos(5  fk)  cos(wk) 

2 

+  45  cos  (apk)  ek  cos(4  fk)  cos(wk)  -  45  ek  cos(4  fk)  cos(wk) 

2  2  2 
+  25  cos  (apk)  ek  cos(3  fk)  cos(wk)  -  25  ek  cos(3  fk)  cos(wk) 

2 

+  60  cos  (apk)  cos(3  fk)  cos(wk)  -  60  cos(3  fk)  cos(wk) 

2 

+  60  cos  (apk)  ek  cos(2  fk)  cos(wk)  -  60  ek  cos(2  fk)  cos(wk) 

2  2  2 
+  30  cos  (apk)  ek  cos(fk)  cos(wk)  -  30  ek  cos(fk)  cos(wk) 

2 

-  60  cos  (apk)  cos(fk)  cos(wk)  +  60  cos(fk)  cos(wk))  zb2/(240  nk) 
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I 


6(ca2S3Sl2)  = 


2  2  2  2 
sin(jk)  n  (24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  24  ek  sin(fk)  cos(4  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  120  ek  sin(fk)  cos(3  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 
+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  104  ek  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  cos  (sk)  sin(wk)  ' 

2  2 

-  160  sin(fk)  cos(2  fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  360  ek  cos(fk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2  2 
+  112  cos  (apk)  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 

-  112  ek  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2  2 
+  80  cos  (apk)  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  80  sin(fk)  cos  (jk)  cos  (wk)  sin(wk) 

2  2 

-  24  cos(apk)  sin(apk)  ek  cos(5  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  120  cos(apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  80  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  160  cos(apk)  sin(apk)  cos(3  fk)  cos(jk)  cos  (wk)  sin(wk) 

2 

-  240  cos(apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2 

-  120  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  cos  (wk)  sin(wk) 

2  2  2 
+  24  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (wk)  sin(wk) 

2  2 
+  120  cos  (apk)  ek  sin(fk)  cos(3  fk)  cos  (wk)  sin(wk) 

2  2  2 
+  104  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2  2 
+  160  cos  (apk)  sin(fk)  cos(2  fk)  cos  (wk)  sin(wk) 

2  2 
+  360  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (wk)  sin(wk) 

2  2  2 
+  112  cos  (apk)  ek  sin(fk)  cos  (wk)  sin(wk) 

2  2 

+  80  cos  (apk)  sin(fk)  cos  (wk)  sin(wk) 

2  2  2 

-  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 
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2  2 
+  6  ek  sin(fk)  cos(4  fk)  cos  (jk)  sin(wk) 

2  2 

-  30  cos  (apk)  ck  sin(fk)  cos(3  fk)  cos  (jk)  sin(yk) 

2 

+  30  ek  sin(fk)  cos(3  fk)  cos  (jk)  sin(wk) 

2  2  2 

-  56  cos  (apk)  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  2 
+  56  ek  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  2 

-  40  cos  (apk)  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2 

+  40  sin(fk)  cos(2  fk)  cos  (jk)  sin(wk) 

2  2 

-  270  cos  (apk)  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2 

+  270  ek  cos(fk)  sin(fk)  cos  (jk)  sin(wk) 

2  2  2 

-  178  cos  (apk)  ek  sin(fk)  cos  (jk)  sin(wk) 

2  2  2  2 
+  178  ek  sin(fk)  cos  (jk)  sin(yk)  -  200  cos  (apk)  sin(fk)  cos  (jk)  sin(wk) 

2  2 
+  200  sin(fk)  cos  (jk)  sin(vk)  +  6  cos(apk)  sin(apk)  ek  cos(5  fk)  cos(jk) 

sin(yk)  +  30  cos(apk)  sin(apk)  ek  cos(4  fk)  cos(jk)  sin(yk) 

2 

+  30  cos(apk)  sin(apk)  ek  cos(3  fk)  cos(jk)  sin(wk) 

+  40  cos(apk)  sin(apk)  cos(3  fk)  cos(jk)  sin(wk) 

+  120  cos(apk)  sin(apk)  ek  cos(2  fk)  cos(jk)  sin(wk) 

2 

+  60  cos(apk)  sin(apk)  ek  cos(fk)  cos(jk)  sin(wk) 

+  120  cos(apk)  sin(apk)  cos(fk)  cos(jk)  sin(wk) 

2  2 

-  6  cos  (apk)  ek  sin(fk)  cos(4  fk)  sin(wk) 

2 

-  30  cos  (apk)  ek  sin(fk)  cos(3  fk)  sin(yk) 

2  2 

-  16  cos  (apk)  ek  sin(fk)  cos(2  fk)  sin(wk) 

2 

-  40  cos  (apk)  sin(fk)  cos(2  fk)  sin(wk) 

2  2  2 

-  30  cos  (apk)  ek  cos(fk)  sin(fk)  sin(wk)  +  22  cos  (apk)  ek  sin(fk)  sin(wk) 

2  2  2  2  3 

+  40  cos  (apk)  sin(fk)  sin(wk)  -  12  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos  (wk) 
2  2  3  2  2 

+  12  ek  cos(5  fk)  cos  (jk)  cos  (wk)  -  60  cos  (apk)  ek  cos(4  fk)  cos  (jk) 

3  2  3 

cos  («k)  +  60  ek  cos(4  fk)  cos  (jk)  cos  (yk) 

2  2  2  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos  (yk) 

2  2  3  2  2 

+  40  ek  cos(3  fk)  cos  (jk)  cos  (yk)  -  80  cos  (apk)  cos(3  fk)  cos  (jk) 

3  2  3 

cos  (yk)  +  80  cos(3  fk)  cos  (jk)  cos  (yk) 
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2  2  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos  (wk) 

2  3  2  2  2 

+  120  ek  cos(2  fk)  cos  (jk)  cos  (wk)  -  60  cos  (apk)  ek  cos(fk)  cos  (jk) 

3  2  2  3 

cos  (wk)  +  60  ek  cos(fk)  cos  (jk)  cos  (wk) 

2  3 

-  48  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos  (wk) 

3 

-  240  cos (apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos  (wk) 

2  3 

-  208  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 

3 

-  320  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos  (wk) 

3 

-  720  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos  (wk) 

2  3 

-  224  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos  (wk) 

3 

-  160  cos(apk)  sin(apk)  sin(fk)  cos(jk)  cos  (wk) 

2  2  3  2  3 

-  12  cos  (apk)  ek  cos(5  fk)  cos  (wk)  -  60  cos  (apk)  ek  cos(4  fk)  cos  (wk) 

2  2  3  2  3 

-  40  cos  (apk)  ek  cos(3  fk)  cos  (wk)  -  80  cos  (apk)  cos(3  fk)  cos  (wk) 

2  3  2  2  3 

-  120  cos  (apk)  ek  cos(2  fk)  cos  (wk)  -  60  cos  (apk)  ek  cos(fk)  cos  (wk) 

2  2  2 
+  9  cos  (apk)  ek  cos(5  fk)  cos  (jk)  cos(wk) 

2  2  2  2 

-  9  ek  cos(5  fk)  cos  (jk)  cos(wk)  +  45  cos  (apk)  ek  cos(4  fk)  cos  (jk) 

2 

cos(wk)  -  45  ek  cos(4  fk)  cos  (jk)  cos(wk) 

2  2  2 

+  45  cos  (apk)  ek  cos(3  fk)  cos  (jk)  cos(wk) 

2  2  2  2 

-  45  ek  cos(3  fk)  cos  (jk)  cos(wk)  +  60  cos  (apk)  cos(3  fk)  cos  (jk)  cos(wk) 

2  2  2 

-  60  cos(3  fk)  cos  (jk)  cos(wk)  +  180  cos  (apk)  ek  cos(2  fk)  cos  (jk)  cos(wk) 

2  2  2  2 

-  180  ek  cos(2  fk)  cos  (jk)  cos(wk)  +  90  cos  (apk)  ek  cos(fk)  cos  (jk) 

2  2 

cos(wk)  *  90  ek  cos(fk)  cos  (jk)  cos(wk) 

2  2  2 
+  180  cos  (apk)  cos(fk)  cos  (jk)  cos(wk)  -  180  cos(fk)  cos  (jk)  cos(wk) 

2 

+  36  cos(apk)  sin(apk)  ek  sin(fk)  cos(4  fk)  cos(jk)  cos(wk) 

+  180  cos(apk)  sin(apk)  ek  sin(fk)  cos(3  fk)  cos(jk)  cos(wk) 

2 

+  176  cos(apk)  sin(apk)  ek  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

+  240  cos(apk)  sin(apk)  sin(fk)  cos(2  fk)  cos(jk)  cos(wk) 

+  660  cos(apk)  sin(apk)  ek  cos(fk)  sin(fk)  cos(jk)  cos(wk) 

2 

+  268  cos(apk)  sin(apk)  ek  sin(fk)  cos(jk)  cos(wk) 

+  240  cos(apk)  sin(apk)  sin(fk)  cos(jk)  cos(wk) 
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2  2  2 
+  9  cob  (apk)  «k  cob(5  Ik)  cot(yk)  +  45  cos  (apk)  «k  cob(4  fk)  coB(wk) 
2  2  2 
+  25  cob  (apk)  «k  cob(3  fk)  coB(wk)  +  60  cob  (apk)  cob(3  fk)  cos(wk) 

2  2  2 
+  60  cob  (apk)  ek  cob(2  fk)  coB(wk)  +  30  cob  (apk)  ek  cos(fk)  cos(irk) 

2 

-  60  cob  (auk)  cos(fk)  coB(wk))  zb2/(240  nk) 
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